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Abstract 


Pionless Effective Field Theory in Few-Nucleon Systems 

A systematic description of low-energy observables in light nuclei is presented. The effective 
held theory formalism without pions is extended to: 

• predictions with next-to-leading-order accuracy for 

the 4-helium binding energy B(a), the triton charge radius, 
and the 3-hehum-neutron scattering length; 

• phase shifts for neutron-deuteron scattering and a-neutron 
low-energy scattering at leading order; 

• the ground states of the 5-helium (with and without Coulomb interaction) 
and 6-hehum isotopes up to next-to-leading order; 

The convergence from leading- to next-to-leading order of the theory is demonstrated for 
correlations between: 

• the triton binding energy B{t) and the triton charge radius; 

• B{t) and the 4-helium binding energy B{a); 

Furthermore, a correlation between B{t) and the scattering length in the singlet S-wave 
channel of neutron-hehum-3 scattering is discovered, and a model-independent estimate 
for the trinucleon binding energy splitting is provided. The results provide evidence for the 
usefulness of the applied power-counting scheme, treating next-to-leading-order interactions 
nonperturbatively and four-nucleon interactions as, at least, one order higher. The 5- 
and 6-hehum ground states are analyzed with a power-counting scheme which includes the 
momentum-dependent next-to-leading order vertices perturbatively. All calculations include 
a full treatment of the Coulomb interaction. The assessment of numerical uncertainties 
associated with the solution of the few-body equation of motion through the Resonating 
Group Method parallels the report of the results for light nuclei in order to establish this 
method as practical for the analysis of systems with up to six particles interacting via 
short-range interactions. 
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[...] 

There is, or was, a poet. His name was Lallafa, and he wrote what are 
widely regarded throughout the Galaxy as being the finest poems in exis¬ 
tence, the “Songs of the Long Land”. 

They are/were unspeakably wonderful. That is to say, you couldn’t speak 
very much of them at once without being so overcome with emotion, truth 
and a sense of wholeness and oneness of things that you wouldn’t pretty 
soon need a brisk walk round the block, possibly pausing at a bar on the 
way back for a quick glass of perspective and soda. They were that good. 
Lallafa had lived in the forests of the Long Lands of Effa. He lived there, 
and he wrote his poems there. He wrote them on pages made of dried habra 
leaves, without the benefit of education or correcting fluid. He wrote about 
the light in the forest and what he thought about that. He wrote about the 
darkness in the forest, and what he thought about that. He wrote about 
the girl who had left him and precisely what he thought about that. 

Long after his death his poems were found and wondered over. News of 
them spread like morning sunlight. For centuries they illuminated and 
watered the lives of many people whose lives might otherwise have been 
darker and drier. 

Then, shortly after the invention of time travel, some major correcting fluid 
manufacturers wondered whether his poems might have been better still if 
he had had access to some high-quality correcting fluid, and whether he 
might be persuaded to say a few words to that effect. 

[...] 


Douglas Adams “Life, the Universe and Everything” 
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I. INTRODUCTION 


Since nucleons are composed of quarks and gluons, the nuclear interaction can be un¬ 
derstood as an effective theory of the strong interaction at a low-energy scale. Effective 
field theories (EFTs), in general, are appropriate representations of particular theories in 
regions of the parameter space of our physical world whereD the latter are assumed to be 
valid but cannot be solved practically. In this sense, Newtonian gravitational dynamics be¬ 
tween two point masses relates General Relativity to a region of small momenta relative to 
the interacting masses. Another example is the classical electrodynamic interaction energy 
between a point charge and a charge distribution - localized within a radius R - which can 
be approximated by the first terms of a multipole expansion for distances ^ R. A quan¬ 
tum EFT, describing properties of electrons (e“) and photons at energies below 1 MeV, is 
Quantum Electrodynamics which does not take into account the weak interaction and its 
accompanying effects of other Leptons. The last example is an EFT with a differing particle 
content, namely e~ and photon gauge fields without fi- and r-leptons and the respective 
neutrinos. The Lagrange density is subject to the same space-time, i.e., Lorentz-, symmetry 
but a 17(1) instead of a SU{2) x U{1) gauge symmetry compared to the underlying weak 
Lagrangean. 

In this work, a “pionless” EFT (EFT(^)EI), with Chiral Perturbation Theory (xPT, see 
[3] for a pedagogical review and references) as the fundamental interaction, is applied to 
few-nucleon systems. This EFT is devised for momenta of the interacting nucleons below 
the pion production threshold and consequently considers solely neutrons and protons as 
degrees of freedom. Similar to the effect of W— and Z—bosons in e~ scattering, pion 
contributions to nucleon amplitudes at low energies can be taken into account by proper 
renormalization. 

The goal of this work is to assess the region of the physical parameter space where 
EFT(if) is applicable. This region is known to include two- and three-nucleon bound- 
and scattering data with and without electromagnetic external probes, but its usefulness 
for the description of amplitudes involving more than three nucleons in the asymptotic 
states is not known. Knowledge about the range of applicability of EFT(^) will further 
the understanding of few-nucleon systems. First, by providing a minimal set of parameters 
for a required accuracy in observables of interest. Therefore, numerical calculations are not 
unnecessarily complicated without gain in accuracy. 


2 


Second, peculiar features of a system, e.g., the unbound five-nucleon state and the halo 
structure of 6-helium, are related to either ^ < 4 parameters (throughout the text, A 
denotes the number of particles in the asymptotic states and B{A) is a binding energy of an 
^-body system which is not necessarily related to a physical state) or to unique properties 
of the respective five- and six-body systems. Analogously, the triton bound state is a 
unique property of the three-nucleon system. Its binding energy, for instance, is in general 
not predicted correctly independently of the accuracy of the description of the two-nucleon 
interaction. In effect, a system of two atoms can behave like a neutron and a proton while 
at the same time the similarity breaks down for observables in the respective three-body 
systems. The three-body properties are unique to the atomic and the nucleon system. Yet 
there are properties in the larger systems which follow from parameters of the interaction 
of the smaller subsystems - little Matryoschka knows her mother to be bigger but cannot 
tell whether she is smiling. We aim to identify universal properties common to systems 
which share certain sets of A < 4 parameters. While universal features of three- and four- 
body systems have been studied theoretically® and experimentally®, universality in A > 5 
systems is new territory. The focus is on bound systems where no a priori arguments about 
applicability or breakdown of the theory are known at the moment. Specifically, no practical 
quantitative procedure is known to properly weigh the high-energy components present in 
the momentum distribution of a bound state, and by that conclude from the distribution of 
nucleon momenta inside the nucleus directly on the breakdown of an EFT. The gravitational 
analog to the nuclei considered here would be a “solar” system of three objects with masses 
and momenta such that bulk properties - average extension and total angular momentum, 
for instance - are to some level described by the Newtonian theory, but the influence of a 
fourth object could potentially lead to an increase of the momenta to an extent that space- 
time curvature can no longer be disregarded, even in the leading approximation. While 
the Newtonian theory is apt for the solar system with its > 8 planets, one may naively 
conclude from the relatively deeply bound four-nucleon state that a description of nuclei 
composed of more than four nucleons in terms of those constituents - namely neutrons 
and protons - alone, i.e., with EFT(^), is inappropriate. Inappropriate in the sense that 
components of a bound state with nucleon momenta larger than the mass of the lightest 
field of the underlying theory that was not included in EFT(|f) - the mass of the pion mn - 
are significant. Significant meaning that the order-by-order removal of parametrically small 
dependencies on renormalization parameters - a feature of all EFTs - cannot be performed. 
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We will provide evidence that this breakdown of the EFT does not occur. 

Related to this puzzle of developing an interaction theory for few-nucleon systems (here, 
few means ^ < 6) systematically from yPT and through that from Quantum Chromo 
Dynamics, is the difficulty of solving the few-body equations of motion. An approach 
refined over decades to yield very accurate predictions of both bound- and scattering nuclear 
observables has been to solve the A—body Schrodinger equation. 


\ i<3 

= |^(D) 


{P^-P3? 

2Amiv 




i<j 


i<j<k 


( 1 ) 


where the center-of-mass motion can be treated independently and the nucleons have the 
same mass m^. The operator Tcm corresponds to the kinetic energy of the center of mass, 
and fj {pj) is the coordinate (momentum) of one of the nucleons in the n-th A-particle 
eigenstate |(;/!)("')). Recent models^® for the potentials Vnn, 3 NI- dots represent possible 4— 
and up to A—body interactions - are based on boson- and meson-exchange phenomenology 
and symmetry considerations. Here, the nonperturbative parts of the EFT(^) interaction 
are related to a potential, too, while perturbative parts are treated as such. The solution to 
eq. Q is obtained numerically, and while one can think of this as a technicality, aloof from 
the theory behind the potential, the method chosen clarifies how uncertainties introduced by 
the numerical technique can be related to the residual renormalization scheme dependence. 

Specifically, our method of choice to solve the Schrodinger equation is a variational 
technique which allows an economic calculation of A < 6 bound and scattering observables. 
The wave function ansatz follows the ideaC®^ of the resonating group structure which 
combines the total wave function of a particular state of a nucleus from components that 
divide the constituent neutrons and protons into separate groups. A hrst approximation 
to the 4-hehum ground state would accordingly be a neutron interacting with a 3-hehum 
core, plus a proton attached to a triton, plus two deuterons bound together. Constraining 
the variational basis with a grouping like this is not done here for A < 4 bound states, for 
which all possible groupings are taken into account - f.e., di-neutron plus di-proton and 
two singlet deuterons in addition to the above mentioned - and a certain grouping merely 
defines the spin- and orbital angular momentum structure. The 6-helium calculations, in 
contrast, will show how a single configuration, namely a 4-helium core with two orbiting 
neutrons, is sufficient for an appropriate description of the ground state. In contrast to 
the most recent incarnation of the method^ which expands the bound fragment states 
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- the triton, 3-helium, and deuteron in the example above - into a shell model basis, the 
method used here employs a nonorthogonal Gaufiian basis allowing for an association of 
single basis vectors with high- or low-energy modes. Therefore, at least a heuristic relation 
of the basis truncation to a regularization of the interaction can be drawn. Loop integrals 
are regularized two times: with a hnite square-integrable basis and the additional vertex 
regulator. Only with an infinite number of width parameters a complete basis would expand 
the momentum-modes responsible for divergences in an amplitude. 

We will utilize the freedom in choosing a regulator for EFT(^) to make the numerical 
calculations highly economical. This efficiency will allow us to analyze correlations amongst 
few-nucleon observables by a variation of a three-body force parameter even for four-nucleon 
scattering observables at next-to-leading order. Correlated observables are a natural con¬ 
sequence of an EFT, which claims to predict low-energy data within defined error margins 
once its coupling constants are htted to other low-energy data. Independent of the specific 
regularization, one set of input data is thus correlated to another, e.g., with an interaction 
fitted to the triton binding energy B{t), predictions for the neutron-deuteron scattering 
length are found in good agreement with data. The approach will be used here for two 
reasons; 

• To reveal potential inconsistencies in the power counting as follows. Assuming a sound 
power counting, a change in the renormalization scheme of the EFT - induced here 
by a variation in the regulator, the input data, and the three-body force - should only 
result in variations of predicted observables consistent with the error estimate of the 
order at which the EFT calculation is conducted. Exceeding this uncertainty margin 
indicates an enhancement of an interaction which was falsely demoted to a higher 
order (see the promotion of a three-nucleon force to leading order when considering 
the correlation between B{t) and '^a^d, see e.g. [T3]i. 

• With a power counting established at a specihc order, discrepancies between (in¬ 
consistent data)/(theoretical predictions) based on different experiments/(interaction 
models or numerical techniques) could be resolved. For example, we will calculate 
the scattering length for elastic 3-hehum-neutron scattering ao(^He-n) correlated to 
B{t). If B{t) is used as fitting input, ao(^He-n) is predicted with an error assessed by 
a change in the renormalization scheme. Any measurement which predicts ao(^He-n) 
outside of this EFT uncertainty margin must be reevaluated. For this method, the 
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knowledge of the “true amount” of the uncertainty is crucial. In the discussion of 
neutron-deuteron scattering we will discuss how various probes of unobservable short- 
distance structure assess uncertainty quite differently. Limiting the probe to a change 
in a regulator parameter, only, will be found to severely overestimate the accuracy of 
EFT(^) at leading order. 

The application of EFT(^) to systems of more than four particles is of particular interest 
because of its implications for other, non-nuclear systems. Peculiar properties associated 
with certain features of the subsystems of nuclei would translate to systems with the same 
number but different species of particles. Specifically, the investigation here aims to resolve 
whether the peculiar unbound hve-nucleon system and the halo structure of the 6-hehum 
bound state are universal properties of systems with unnaturally large two-body scattering 
lengths and a three-body bound state. Such a finding could be related to atomic systems, 
where experimentally^ similar conditions regarding as^t and B{t) can be realized. 

In this work, we study nuclear systems only, while the results can be related to other 
systems which interact via forces of the same universality class. I start with a general 
overview of the EFT formalism. 

II. NUCLEAR EFFECTIVE FIELD THEORY 

Generalities. An effective field theory exploits the existence of different scales 
in a system. While an underlying theory is valid at least up to an energy scale M with 
a Lagrange density T ('k) constrained by spacetime and internal symmetries, the EFT is 
useful only for momenta Q <C M. In general, either the ignorance about the underlying 
theory, or the inability/impracticality to solve it in the energy interval of interest, motivate 
the model-independent EFT approach. The partition function of the effective theory, 

Z = , (2) 

from which the S'-matrix is derived, depends on a different set of degrees of freedom (DoF) 
tp, namely either only low-momentum modes of T or composed states of the fundamental 
fields, for example, nucleons instead of quarks or a particles instead of pairs of neutrons and 
protons. The symmetry group is the same or an approximation of the original, constraining 
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the operators Q with mass dimension <5 in the effective Lagrangean, 

N{5) 

£('») t'') a I'*) • (3) 

S<Sm i{S) 

If expressed in the same DoF, the S'-matrix of the underlying theory is recovered for dm —t oo 
and momenta of order Q. The low-energy coupling constants (LEG) gi depend on a renor¬ 
malization parameter A which specifies a regularization scheme and renormalization condi¬ 
tions and thus the balance of loops and vertices for a specific contribution to the S'-matrix. 
Assuming a natural size for the LECs (see eq. ([^), a truncation of the sum in eq. Q at 
A^max will yield a Lagrange density which contributes all terms, up to a certain power, of 
an expansion of the S-matrix in the small parameter Q/M. In addition to this ordering 
scheme of the interaction terms in £ - the vertex expansion - the EET framework requires 
an ordering of iterations of those operators - the loop expansion. 

Two examples of phenomena determining the scale M are: 

• the mass of a particle associated with a DoF of the EFT whose mass is too large 
relative to its momentum to be produced or to decay; this is relevant in deriving the 
kinetic-energy term for nucleons below; 

• a scalar field operator composed of DoF from the underlying the¬ 
ory with a non-zero vacuum expectation value; relevant here: 

T G {neutron(n), proton(p),pion(7r), A — resonance}. 

Having set M, the theory is expected to be useful for observables which are dominated by 
external momenta Q <C M, and to break down if Q is of the same magnitude or larger than 
M. The contributions from the most general Lagrangean compatible with the symmetries 
of the system, conjectured to yield the most general 5-matrix with the same symmetries for 
Q <C M, are subjected to an ordering scheme in powers of Q. This scheme is to be verified 
a posteriori by the analysis of the dependence of low-energy observables on the renormal¬ 
ization scheme parameterized by A. If this dependence exceeds the contributions expected 
from operators with N > Am ax . the power counting was based on false assumptions, and 
one has to consider, e.g.: 

• LEC(s) of unnatural size; 

• Typical momenta in a system which are too large to put a certain observable within 
the range of applicability of the EFT. 
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TABLE I: Data for light nuclei. To estimate the low-momentum scale Ton = 938.858 MeV is used. 



deuteron(d) 

triton(^H,t) 

3-lielium(®He) 4-helium(a) 6-helmm(®He) 

B 

[MeV] 

2.2245727(20)Ii^ 

8.481855(13)Ii^ 

7.718109(10)Ii^ 28.295875(26) 29.2712(10)1^^ 

{rXC 

[fm] 

2.0952(60)^1 

1.63(3)IiS] 

1.976(15)11^ 1.671(14)11^ 2.054(14)123 


[MeV] 

45.70 

72.86 

69.50 115.3 117.2 


Degrees of Freedom. The objects of interest are nuclear systems with A < 6. 
For scattering observables the interest is in center of mass energies Ecm < 10 MeV. The 
naive estimates for Q in the bound systems considered in this work together with relevant 
data are shown in table |Ij Nuclear processes at momenta of order Q do not probe details of 
the internal structure of the nucleons, i.e., quarks and gluons. Transitions to excited states 
such as the A or the decay into lighter particles, e.g., vr- or p-mesons, with high momenta, 
are strongly suppressed. With Q < (mass of the pion), even virtual pion exchanges can 
be approximated by contact interactions. As a consequence, the theory is formulated with 
non-relativistic fermions, namely neutrons and protons, as the sole DoF. The number of 
nucleons A in the incoming state is conserved and only affected by interactions represented 
by operators with < 2A nucleon fields. The scales of the “pionless” effective field theory 
(EFT(^)EI]) are: 

• high-momentum-/breakdown-scale: M ~ tUtt 

• low-momentum scale: Q 

The low-momentum scale is estimated here by the assumed typical momentum of a nucleon 
in the deuteron with binding energy B[d). For scattering calculations, the effective theory 
is useful for external momenta up to this order. In bound systems, the usefulness depends 
on the momentum distribution of individual nucleons. In general, the bound state is a 
superposition of states with relative nucleon momenta ranging from zero to infinity. Only 
if the momentum distribution function supresses states corresponding to momenta larger 
than Q sufficiently, a useful description of the bound system can be expected. For the 
deuteron, the formula in the last line of table estimates typical momenta of the nucleons 
in it well, if the high-momentum states are strongly suppressed. For larger systems the 
angular momentum structure is more complicated. It is, for instance, possible to find 









two of four particles in a relative D-wave state while the other two are in a relative S- 
wave. Hence, the naive estimate for a typical momentum given in table |I] is no a priori 
argument for the breakdown of EFT(y') in the respective system. For a defined external 
momentum, however, the expansion of observables provided by this theory is appropriate 
for Q <C C>(m 7 r), becomes less useful for Q < 0{mT^), and is not applicable for Q > 0 ( 171 ^^). 

Symmetries. Having defined the DoF, the interaction terms in the Lagrange den¬ 
sity are constrained by spacetime and internal symmetries. Nucleons with small momenta 
relative to their rest mass can be approximated by Pauli instead of Dirac spinors, and the 
Lagrangean should only be invariant under Galilei transformations and not under elements 
of the proper Lorentz group. Parity remains as a symmetry of this theory for the strong 
interaction. The transition from a relativistic formulation of the kinetic energy density to 
its low-energy version is summarized following®!^, to make the connection between the 
relativistic and non-relativistic theory. 

{i0 - mj,)n ^ £kin = ^ = '^Nl'j/!v^d''Ny 

V V 

(4) 

1. Fermion fields n(neutron or proton) with 4-momentum pi, = m^Vi, do not change 

their velocity in the limit —)■ 00 under a momentum transfer which is small 

relative to m^. Hence, modes with different Vy decouple. 

2. The decoupling makes a velocity-dependent field transformation 

convenient. The fields are split into particle annihilation and anti-particle creation 
operators -t- N~. The momentum transfer is assumed to be small relative to 

the rest mass, not providing enough energy for the creation of particle-anti-particle 
pairs: $N~ = = 0. 

3. For small momenta the Lagrangean must be invariant under a field transformation 

resembling a small boost, Ny —Ny = (only here, N" is a scalar field). 

The effective Lagrangean up to O ( — ) is given. 
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4. The final form of the kinetic term of the effective Lagrange density, used here, is 
obtained in the rest frame of the nucleon after field transformations involving ivyd^ 


and 9^. The operator N = 


P 


n 


holds the neutrons and protons chosen each as a 


two component Pauli spinor with mass dimension [A^] = |. The corrections from the 


mass difference are assumed to be less than order m 


-2 


The internal symmetry transformations act on the two components of the isovector 

(p\ 

= I I . In this work, the only interaction which breaks the isospin symmetry - i.e., 
distinguishes between two neutrons and two protons - is the Coulomb term. 

Power Counting. The EFT formalism orders contributions to an observable in 
powers of the typical external momentum Q. Initially, the ordering scheme is based on the 
assumption that all coupling constants are of natural size, i.e.: 


Ci (^) 

gi (A) = - 4 : with dimensionless Ci(l) of order 1 . 


(5) 


This implies an ordering scheme with a truncated expansion £('5<^max) (vertex expansion), 
and also a limited number of iterations (loop expansion) of the interaction terms contained 
in C. Operators of higher mass dimension or higher iterations of terms with S < dmax would 

/ Q \ ^(^max) 

be of order ( ^ j and therefore suppressed. For dmax = 8 the effective Lagrange 

density is: 


£(8) ^ +- JV + C}^^{N^N){N^N) + C^^{N^a,N){N^aiN) 

\ 2mN J 

-IcNLo [(7vta^7v)2 + 

- (N'fd.N) (d^N'<)N + [N'idfN + dfN'iN] 

* xxNLO, 


[(iVta^iv)[(a,iVt)afciV] + [(d,lVt)iV](iVV,dfe7V)] 

-(iv^iv) [(a,ivt)a,dfeivl + (ivV.Tv) [(a,ivt)dfe7v^ ^ ^ ^ 


41 


NLO I 


+ -C. 


4-r 




+ ( 5,,5ki 


[{d,djN^)akN] + {N^akd.d.N)] (N^aiN) 


2 1 
1 /I 


iS.kSji + 6u6kj) + iN^akd,N)[{d,N^)aiN] 


~8 ~ " 3 CrLO)% 4 i ) [{d^N^akd.N) 

+ idjN^akm)]iN^aiN) . 


( 6 ) 
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TABLE II; Experimentally neutron-proton effective range parameter for the spin singlet (s) and 
triplet (t) S-wave channel. The conversion factor is fic = 1 = 197.316132894 MeV-fm. 


Os at rs rt 


[fm] -23.748(10) 5.424(4) 2.75(5) 1.759(5) 
[MeV-i] -0.12036 0.02749 0.0139 0.008915 


The vector a consists of the Pauli spin matrices, lower case Arabic letters label vector 
components in a Cartesian basis, and the sum convention is understood. In leading order 
(LO), only operators with d < 6 - namely C\^ - are included, while next-to-leading order 
(NLO) means also to consider 5 = 8. To retain the S'-matrix for low-energy two-nucleon 
reactions, the assnmption that all LECs are of natnral size has to be given up, 

because the resulting perturbative treatment of all interactions cannot reproduce poles of 
the amplitude at momenta Q <C M. The poles, if there are any, of an amplitude derived 
from a theory with natural coupling constants only, are all at momenta Q > M. The 
peculiar analytic structure of the amplitude is resembled in the relative size of the effective 
rangel^ parameters in both the singlet- and triplet-neutron-proton (total spin s = 0,1) 
channel (table 0 The scattering lengths Us^t are considerably larger than the effective 
ranges which are 0{1/M). Hence, a<j and at constitute additional low-momentum 
scales, associated with a virtual and a real bound state close to the scattering threshold. 

^ <C M ~ O(m^) . (7) 

I I 

With the additional scale a coupling constant of size < 71^2 (A) = with 

Cl, 2 (1) rsj 0(1) does no longer justify a truncation of the loop expansion. This assump¬ 
tion is made for C\^. Because of this unnatural size, an iteration of such an interaction 
does not change the order of the respective contribution. Diagrammatically, this is shown 
in fig. where all summands in the first row are of the same order with only momentum- 
independent vertices, while insertions of vertices with natnral sizes decrease the size of the 
contribution to the amplitude. 
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loop 



FIG. 1: Lowest order contributions to the amplitude T. The vertices are labeled by the mass 
dimension of the operator. The first line of diagrams is 


The leading-order amplitude can be calculated analyticallyESl in momentum space: 

d^l dlo i i 


T{p,p',E) = -iC + i-iCf 
= -iC + {-iC)^ 
= -iC + i-iC)"^ 


(27r)3 Itt 1 , TP P I / P I 

V J lo + E - ^ -b ^e -Lo - ^ + ^e 


+ o{c^) 


dH 1 


(27ri) 


-1 


(27r)3 27r^ pj _^ _l_ 


mN 


d^t 


— I 


(27r)3^_ A + ^ 


-+ o(A) 

ie 

+ o(A) 


mN 


le 


=-j/A = (-27r)-3 r^dl - 

' / JO E-l^+iE 




1 


A - — 


n=0 




( 8 ) 


The constant C is a placeholder for the propagator is derived from the kinetic part 

-Q_ 

of £, the total kinetic energy in the center of mass system \s E = and the subscript A 
specihes a regularization scheme for the divergent integral. 
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Parameter Determination. Changes in the high-energy/short-distance struc¬ 
ture of the theory, performed either by integrating out DoF or by the regularization, man¬ 
ifest themself in the values of the LECs. At the moment, a direct calculation of the LECs 
from the underlying theory is not feasible. Instead, low-energy data related to coefficients of 
a power expansion of the amplitude is used to determine the LECs. The expansion sketched 
in fig. [^is identical to the effective range expansion ^3. The latter approximates the scat¬ 
tering amplitude (parameterized by a phase shift 6 i), for a relative angular momentum I 
and center of mass momentum p = 


I 


coi 5i = -h 

ai 2 


+ 


OO 

E (?^) 2n 

v\ ’p 

71=2 


0 (S-wave) : T{p) = 


dvr 1 
ttin p cot 6 o — ip 


(9) 


This expansion is around zero momentum and introduces the coefficients ai (scattering 
length), ri (effective range), and vi. Before relating this expression to the amplitude from 
eq. ([^, a regularization scheme has to be specified. With the momentum dependence of I\ 
thus fixed, the result in eq. Q is expanded in powers of p to obtain a form analogous to 
eq. ([^. Equation]^ is of Lippman-Schwinger-type with a potential V{p,p') = —iC. With a 
regulator function 


f lim f{p, A) = 1 
f{p, A) such that < 


I /(P,A) 


( 10 ) 


1 for \p\ of order Q and A > 0{M) 


and using the separable form 


Vr{p,p',K) = f{p,A){-iC)f{P,A) 


( 11 ) 


in eq. pj), yields: 




-i-iA 


Ja = dvr 


dl Pp{l,A) 


( 12 ) 


To have a GauBian regulator in coordinate space, we chose a different nonseparable regulator 


which will be explained below (eq. 20). A power expansion in \p\ of the integral I a can 
now be related to the effective range amplitude. In leading order, EFT(|f) reproduces the 
experimental scattering lengths (singlet: 5 = 0, triplet: 5 = 1) for 

dTT 


C\^ + (d5 - 3) = 


LO 


mNvr\/2 — rrijqAy/^as^t 


(13) 
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Having determined the LECs Ci^, the corresponding potential can be used in the 
Schrodinger equation - equivalent to the Lippman Schwinger equation - to make predic¬ 
tions. At this point, I summarize the prescription which starts from a Lagrange density C 
and results in a potential that is used in a multi-body Schrodinger equation. 

1. The first part of the power-counting scheme is a truncation of the vertex expansion 
of C based on dimensional analysis and natural coupling constants. I will consider 
two orders in this work: 

• leading order (LO): £('5max=6) 

• next-to-leading order (NLO): £('5max=8) 

2. At each order, the vertex structure of C is identified with a potential, which is given 
in momentum space: 

LO: V^§{q,k)= C\^ + Cl^ai-a2 ; (14) 

NLO : k) = + C^^^k^ + <ti • 0^2 (15) 

. g X fc + 0^0^-. ^ . 

Where q = p — p' and k = 2^ is used. 

3. The second part of the power-counting scheme is an estimate of the relative size of 
terms in the loop expansion. This estimate determines whether or not an operator 
in the potential can be treated as a perturbation. 

4. The theory’s renormalization scheme is defined by a regulator function and the low- 
energy data used to fit the LECs. 


5. Predictions of the resulting interaction are expected to be consistent with data within 


the uncertainty margins: 

A(LO) ~ O fPS ^ PS ^ , (16) 

^ ^ \M J rUn 3 ' ^ ^ 

A(NLO) ~ O (g) » i . 


where a typical momentum of a nucleon within the deuteron of 7 t = 
is assumed. Larger deviations resemble either an insufficient power counting, e.g., 
LEC(s) of unnatural size, or a failure of the EFT because the typical low-energy scale 
of the system is larger than expected. 
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Renormalization Scheme Dependence. I will use two related criteria to vali¬ 
date the power-counting scheme for a specific observable: 

• order-by-order convergence of the prediction consistent with an expansion parameter 



• renormalization scheme dependence of < O (|) at LO and < O (|) at NLO, i.e., an 


maximal deviation of approximately 30% (10%) from the experimental value at LO 
(NLO). 


In addition to the regulator value A and a subset T)' <Z T) oi all low-energy “P(ata)”, I 


define the cutoff parameter Arqm (see ch. IV) as a third parameter to completely specify 


the renormalization scheme. Renormalization, in general, constitutes a modihcation of 


unobservable modes, i.e., with momenta outside the range of validity of the theory, while 


imposing low-energy constraints - renormalization conditions. The lower bound for A was 
chosen of the order of the breakdown scale, Amin = 150 MeV, in order not to affect the 
interaction for center-of-mass momenta less than 100 MeV. Argm parameterizes a change of 
high-energy modes by means of adding or excluding GauBian basis functions with support for 
|f| <C For given A, Argm, and V’ , a correct power counting predicts o €V\V’ with an 
uncertainty depending on the order of the calculation. A different renormalization scheme 
will affect those predictions only within this uncertainty range. A larger deviation, i.e., 
larger scheme dependence, indicates an insufficient power-counting scheme. A prominent 
example is the promotion of a three-body interaction to leading order in the triton channel 
(s = |, with the S-wave neutron-deuteron scattering length denoted by a(nd)t)- Considering 
V = {a^, at, B{t), a(nd)t} and V’ = {a^, at] to fit the predictions for T>\V’ will show 

a strong dependence on A, known as the Phillips lineE3 . Promoting the <5 = 9 operator. 


£^=9 = (Art N) (Vt TiN) (Vt TiN) , 


(17) 


to leading order, and using V = yields predictions for a(„rf)t consistent with 

the LO uncertainty estimate, while not affecting A = 2 observables El EH ESI _ Accordingly, 


I extend the leading-order potential in eq. (14) for calculations in A > 3 systems by the 


operator: 


^Mi(P12,P23,P'i2jP'23) — C'sNi'R ‘ ^2 


(18) 
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Universality. The unnatural scaling of Cgj^j is a consequence of a finite quan¬ 
tity characterizing the three-body system in the limit of a scale-free two-body subsystem 
in which, e.g., diverging two-fermion scattering lengths, ag^t oo. This property was 
discovered®^ as the universal ratio of binding energies of successive states in the spectrum 
of three identical bosons. In the limit |a| —)• oo, this spectrum consists of an infinite number 
of bound states with energies which accumulate at the lowest break-up threshold. 

These binding energies approach a universal ratio given by: 

^(n+l)(3) 


lim , ,, , 

n->-oo i?(»^)(3) 


515.03 


for |a| —oo 


(19) 


Here, n = 1 labels the largest binding energy, and universal refers to the independence 
of that ratio with regards to the specihc form of the two-body interaction that generates 


the diverging scattering length. For instance, the LECs in eq. (14) can be tuned such that 
—)• oo for different values of A, i.e., different shapes of the corresponding potentials. 


However, the ratio given in eq. (19) will be independent of A, it is universal. In contrast, the 
absolute position of the spectrum, e.g., fixed by the smallest eigenvalue in the three-body 
spectrum, is no universal consequence of the infinite scattering length - it differs in general 
for different A. If, however, this eigenvalue is fixed, so is the complete spectrum, which is 
why in addition to only one additional parameter is needed to select the system that 
corresponds to the triton from the universality class, associated with In practice, two 
approximations result in a finite three-body spectrum: First, the scattering lengths are 
finite, yet large relative to a range set by the breakdown scale of the EFT, and second, the 
spectrum is calculated in a finite basis. 

With the extended set V = {ag,at, the predictions of the EFT for V\V' are 

expected to by accurate up to corrections of the order of about 30 %. Strong evidence was 
found recentlyED that this assertion holds even if T> includes low-energy four-body data. 
In contrast to the three-nucleon system, the four-nucleon system is then at this order a 
universal consequence of its constituent subsystems without universal scales. An analogous 
promotion of a four-body counter term to LO is therefore not necessary. 
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III. NUCLEAR POTENTIALS FROM EFT 


I will derive three classes of nucleon-nucleon (NN) potentials. The first class consists 
of the leading-order interactions, while the second and third class are both of next-to- 
leading order accuracy but include different subsets of higher-order contributions. With 
the potentials, A-body observables are predicted from solutions of the Schrodinger equation 
eq. Q. The eigenstate is obtained with a variational method approximating the wave 
function in coordinate space. This method is based on the idea of resonantly rearranging 
groups of particles within the A-body system and is reviewed together with the extensions 
for its application to the EFT(^) potentials in ch. |IV[ The necessary Fourier transform of 
the NLO operators in eq. © is given in app. but the form of the regulator shall be 
discussed here with the LO potential as an example. To minimize the numerical effort, I 
choose a momentum space regulator that transforms into a Gaufiian function on the relative 
coordinate between the interacting particles. The conjugate momentum to this coordinate 
r (r') before (after) the interaction is p {p'), fig. Choosing the same regulator for all 
vertices, 

fip,p',A) = U(p,p',A) =, (20) 


yields the convenient leading-order form: 
1 




(27r)6 

1 

(27r)6 

1 


d^p'd^p {r'\p'){p'\V\p){p\r) = 


1 

(27r) 

d^kSq V(fc, 


d^p'd^p ,p)e-^^-P 


d?kd^q 

3 4'' 2 ' n( r — ' 


(27r)® j 




Stts 


( 21 ) 


A“5((r^-hr0/2) 

The momentum-space regulator is hence not separable in p and p', a property used to derive 


the solution of eq. (12). Consequently, this solution provides only a useful starting value 
for the numerical fit of the LO LFCs, as will be shown below. 
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Leading-Order Potential. At leading order, the theory is based, in the two- 
nucleon sector, on Iterations of the two vertices are all of the same order due to the 
additional unnatural scales of order Q in both NN S-wave channels. This power counting 
is therefore correctly implemented by solving the Schrodinger equation with a potential 

« = . (22) 

i<j V87r2 / 

The two independent LECs, are fitted to the scattering length in the spin-0 S-wave 
neutron-proton (np) channel and to either at(ripiet) or the deuteron binding energy B{d) in 
the spin-1 S-wave np channel. B{d) is identified with a pole of the amplitude, which is in 
LO of the effective range expansion: 




2 

t 


m-^a^ 


(23) 


B[d)\ 

7 ^lat=. 


dt —flt,exp 


1.41 MeV or ad 


B{d)=B(d,exp) 


0.0219 MeV 


-1 


Using the experimental deuteron binding energy B{d, exp) as input yields a scattering length 
that deviates from data within the LO uncertainty range. B{d) is consequently identified 


with a low-energy observable that can also be used to determine LO LECs. Table III 


lists the leading-order potentials derived for this work, where B(t) and B(a) were calcu¬ 
lated with = 0. Eor A < 800 MeV, the dependence of the binding energies on the 
ROM model space was found small: < q g]^ (models spaces wl2(0) are 

defined in table VII & VIII &: IX), while for A = 1600 MeV only wl2{) is to be used 
because of the shorter range of this interaction. In practice, the scattering lengths are 
calculated from the phase shift 5 in the respective channel at an energy Eq = 0.4 keV via 

two-nucleon sector, the phase shifts were calculated by numer¬ 
ically integrating the Schrodinger equation using Mathematica and alternatively with the 
Resonating Group Method (RGM). The former uses a steepest descent optimization method 
to fit the LEG, while the latter employs a genetic algorithm (app. 0. The value predicted 


by eq. (13) for the separable regulator is used either as a starting value (steepest descent) or 
center of the search interval (genetic algorithm) in the two approaches. The LEGs obtained 
with both methods differ by <C 1%. This difference transients to a numerical uncertainty 
for predictions of observables which is insignificant relative to the expected systematic error 
at leading order. 
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TABLE III: Parameters of leading-order potentials of the form eq. (211. If unspecified, quantities 
are given in MeV. 


A [MeV] 



8"V-iA®C2° 

[Csl \Ct\ 

A® 

V 

B{d) 

Bit) 

Bia) 

200 

0.25682 

-41.549 

-5.1372 

0.31 

2.6 

0.56 


2.5315 

9.1529 

31.975 


0.25685 

-40.592 

-4.8183 

0.38 

2.5 

0.028 

as, B{d) 

2.2245 

8.5156 

30.371 

300 

0.57784 

-80.727 

-6.7272 

0.55 

1.4 

5.83 


2.0587 

13.047 

53.308 


0.57791 

-81.480 

-6.9786 

0.54 

1.5 

6.47 

as, B{d) 

2.2246 

13.476 

54.459 

800 

4.1091 

-493.38 

-15.467 

0.34 

0.44 

190.0 

0>s,t 

1.6168 

54.285 

262.59 


4.1096 

-500.34 

-17.789 

0.33 

0.45 

221.5 

as, B{d) 

2.2246 

57.539 

271.47 

1600 

16.436 

-1897.8 

-29.731 

0.18 

0.21 

2792.0 


1.5145 

187.288 

374.508 


16.438 

-1913.8 

-35.091 

0.18 

0.21 

4551.0 

as, B(d) 

2.2246 

194.24 

382.603 


The leading-order LECs we found are of a size consistent with a basic EFT assumption 
(discussion of eq. §)• Specifically, oneE3 estimates, 


+ (45 - 3)C^ 


LO 


dvr 


7 

^O(l) 


10“^ MeV“2 


( 24 ) 


where I set M and the exact value of 7 depends on the chosen renormalization 

scheme, e.g., parameterized by A. The value of eq. [^is not exceeded by the numbers for 
Cs[T) = 10 ‘^- 87 r 5 A“^(C }"‘^-|-(45 — 3 )C' 2 °) as given in table III which are consitently smaller 
than 10“'^ MeV“^. 

The two-body interaction thus determined, without proper renormalization through a 
three-body vertex, predicts a triton binding energy B{t) that increases with A, i.e., with 
decreasing range of the potential at constant B[d). This behavior is known as the Thomas 
effect E2]. This divergence disappears naturally with an additional renormalization con¬ 
straint at LO on a three-body observable. I choose a three-body interaction of the form 
in eq. ( 0 : regularized analogous to the two-nucleon interaction to yield an operator 
coordinate representation of the form: 

A 


m 


^LO _ 

'^3NI — 2^ 


i<j<k 

cyclic 


A3 

3 ' 

87r2 


AifS 
4 ^j 


A3 

3 

87r2 




^7^ a 


''3NI Ai • Tj 


( 25 ) 
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In addition to the the physical constraint put on the minimal value for A, practical cal¬ 
culations were facilitated by also imposing an upper bound. The increasing B{t) requires 
a correspondingly stronger 3NI to fit data and poses an inconvenience regarding the nu¬ 


merical technique (see ch. IV for terminology). Specifically, the variational basis has to 
expand two wave functions which differ widely in their spacial extension, with the more 
deeply bound system naturally confined to a less extended region in space (see discussion 


of the triton charge radius in ch. VIA). Hence, the basis has to expand a wave function 
with small support - corresponding to the LO state without 3NI - relative to the support 
of the renormalized state. By using cutoff values in the range 300 MeV < A < 800 MeV 
for which B{t) at LO without 3NI does not exceed 60 MeV, a single triton basis could be 
used for all cutoffs. This interval does not present a principal limitation on the regulator 
imposed by the ROM. For As exceeding the boundaries, the basis must be augmented with 
sufficiently narrow width parameters to account for the structure of the wave function at 
distances r ~ A“^. Albeit constraining the cutoff to this interval, the reduction of the basis 
size for a triton substructure within A > 3 trial wave functions has to be performed for 
each A separately (see discussion of the ^He core as a subsystem of ®He in ch. 
analogous case). 


IX 


for an 


An elegant way of dealing with larger As is provided by the fact that the number of 
bound states increases with the cutoff. This opens the possibility to ht not the deep ground 
state to the triton but rather an excited state - which is already closer to the datum - 
. The advantage is a relatively small contribution from the 3NI at the expense of an 
assessment of the convergence of the excited state within the RGM model space. Although 
this convergence is guaranteed (see Mini-Max theorem, e.g., [M])) in practice, a larger basis 
are necessary. For this reason I do not draw this option and choose A from the above defined 
interval, where for all practical bases only one bound three-nucleon state exists. In ch. |IX[ 
this identification of a physical nucleon with an excited state in the spectrum of EFT(^) 
is discussed further. There, in the six-nucleon system, we will find multiple bound states. 
If any of them corresponds to the physical 6-helium ground state or if the the theory does 
not sustain a shallow bound state at all at the respective order will be a central problem 
in that analysis. Here, a strong cutoff dependence observed in NLO corrections for those 
deep states is suggested as only one discriminating factor between states in- and outside 
the range of validity of EFT(^). 
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TABLE IV: Next-to-leading order constraints on observables receiving dominant contributions from 
terms of higher order in the vertex expansion. 



nonpert. 

pert. 


C(0.1-,5fWA) 

0 

ei ((35 i|VNLO|3^^^^ 

C (0.1 • efWA) 

0 

Next-to-Leading-Order Potential. 

Truncating 

the 


<^max = 8 with a promoted three-nucleon interaction (3NI), the numerical uncertainty in 
the amplitude - calculated as sketched in hg. [T]- is estimated to decrease by a relative factor 
of The effective range expansion provides two additional parameters of 0{M~^). Specif¬ 
ically, the NLO input data is taken as a 4-element subset oiV = {og, ot, S(t)}. 

The tree-level amplitude of is given by eq. (15), and transforms (app. [a| with a regu¬ 


lator function of the form eq. (20) into a coordinate space function: 

A 

yNLO _ /q (A, r) {Ai + ■ (Jj) + (A3 -I- A4^(Ti ■ (jj) |/o (A, r), V^| 


i<j 

-t/o (A, r) (As -I- A^Si ■ aj) + Iq (A, r) A7L ■ S + Iq (A, r) Ag ( • mj • f- 


-Ag -;/q (A,r) , [S’* (g) 5®] ^ (g) [erf (g) cr|] ‘ 


n 0 


(26) 


-» 1 ^ ^ -2 
with L = —irx\7 , S —-{ai + aj) , Ig (A, r) = e“~’' 


and nine redefined LECs Ai which are functions of the original xhe explicit form 


of this dependence on jg given in table XI in the appendix. Contributions from 


EFT (it) to effective range parameters for I = 1 (P-wave) and the transition / = 0 —t t = 2 
(SD-mixing) are suppressed relative to contributions to the effective ranges rg^t by at least 
another factor of This allows for an arbitrary constraint on the predictions of 
in these partial waves, because only the next orders will introduce LECs that have to be 
fitted to physical values. I will use two types of constraints, depending on whether the 
NLO operators are iterated (nonperturbative treatment) or not. The specific constraints 


are given in table IV We chose the values for the phase shifts 61 and ei as small fractions 
of the Nijmegen partial wave data^ (PWA) because then LECs are of magnitudes that do 
not pose numerical problems. 
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Nonperturbative (nNLO). 


The first scheme I used, identifies the functions given in eq. (26) & (25) with a potential, 
from which predictions are made by solving the corresponding Schrodinger equation. This 
equation of motion for the wave function can be transformed into a Lippman Schwinger 
type integral equation like eq. ([^, and by that it can be seen that its solution is found by 
iterating all operators in the potential. The amplitude, thus obtained, includes in addition 
to the LO and NLO diagrams shown in fig. contributions of O The effect is 

analogous to that of the regulator function. The resulting amplitude can be cast into the 
form 


nipi )« 


TJCT + (<^1 + fi) + *P 


(27) 




i=0 


with expansion coefficients di ^ 0 even for i > 1 resulting from the iterated NLO vertices. 


The fi’s parameterize the expansion of the regularization analogous to I a in eq. (12) at 
LO. The cutoff dependence leaking through non-zero coefficients with i > 1 has to be 
< O for this power-counting scheme to be useful. The di will be a linear combination 


of LECs from eq. (15), and therefore, this criterion could be met by a cancellation of 
large coefficients. Although this is acceptable from the EFT point of view, the associated 
difference in the relative magnitude of matrix elements causes numerical uncertainties which 
increase significantly beyond A = 4. For A < 4, however, explicit calculations will provide 
evidence that the criterion is met, and the numerical treatment is consistent in this respect 
with calculations using different renormalization schemes for which analytical solutions exist 
even at NLO. 

The redefined LECs are chosen with the genetic algorithm to minimize the weight func¬ 
tion 




E 


o(exp) — o(trial)| 
|o(exp)| 


(28) 


with P' = { (Re,Im{e2™(^)},TO e {d (I’^S'o.i,^’^ Po. 1 . 2 , ei)}, E g]0 MeV,0.35 MeV]) ,E(d)| 

In order to assess convergence and accuracy, potentials with different short distance struc¬ 
ture - parameterized by the regulator cutoff A and the input data - are derived (table |V]) . 
The predicted NN phase shifts, shown in fig. are consistent with expectations regarding 
LO to NLO convergence and accuracy at the respective order. In particular, the validity of 
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FIG. 2: Selected neutron-proton scattering phase shifts for the potential and values from the 
Nijmegen partial wave analysis^S which were used to fit the LECs. 


the constraint on the P-wave phase shifts is apparently satished and the predictions for en¬ 
ergies larger than the fitting interval remain small. We iterate that those constraints reduce 


the number of independent LECs in eq. (26) to five. Also, as argued above, they are consis¬ 
tent with the EFT philosophy that higher order interactions can only induce higher order 
corrections in observables. Different P-wave constraints, corresponding to either option in 


table IV, allow for such a controlled modification of short-distance structure. Enforcing zero 
P-wave phase shifts provides no significant gain in computer time in calculations in A < 4 
systems. We therefore exploit the additional handle on higher order effects provided by the 
first option which allows non-zero values to gauge the accuracy of a NLO calculation. A 
consequence of this approach are non-vanishing rank one and two interactions correspond¬ 
ing to LECs Ay and Ag^g, respectively. These tensor structures are also found when one 
implements the SD-interactions, which enter beyond NLO. We choose to include these as 
representations of higher-order effects and to constrain their parameters by the SD-mixing 
angle. 
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TABLE V: Numerical values of the cutoff A and the LECs of the NLO EFT(y-) potentials. The 
LECs were fitted to B{d,exp) [3S] and NN phase shifts. 


A 

[MeV] 

Ai 

[MeV] 

A 2 

[MeV] 

As 

[MeV-fm^j 

A4 

[MeV-fm^] 

As 

[MeV.fm“^] 

Aq 

[MeV-fm~^] 

A7 

[MeV] 

As 

[MeV.fm“^] 

Ag 

[MeV-fm^] 


414 

-143 

464 

-57.0 

68.8 

-52.4 

-346 

-3.97 

-62.1 

-0.111 


432 

-612 

944 

-157 

110 

351 

-723 

66.0 

-168 

-0.137 


544 

-1224 

1036 

-432 

336 

1704 

-1851 

-78.9 

-378 

-0.130 

% 

544 

-1625 

-89.6 

-298 

49.8 

1870 

34.3 

-990 

-734 

2.13 • 10“® 

% 

544 

-405 

317 

-125 

108 

287 

-493 

-271 

-360 

-9.88- 10“'‘ 


648 

84.9 

-324 

-888 

404 

4845 

-2064 

-342 

-2319 

-0.716 


648 

-1316 

1039 

-530 

431 

3775 

-3666 

-75.4 

-1052 

-0.178 

00 

672 

301 

-143 

-708 

480 

3527 

-3450 

-192 

-1613 

-0.257 


672 

-158 

-201 

-362 

275 

2001 

-1659 

79.7 

-1270 

-0.0850 


Perturbative (pNLO). 

Although the partial inclusion of higher-order contributions through an iteration of NLO 
operators in the course of solving the Schrodinger equation does not result in a cutoff depen¬ 
dence inconsistent with the expected NLO uncertainty, i.e., the uncertainty and convergence 
estimates are still valid, a perturbative treatment is expected to have two main advantages: 


1. The numerical problems related to the cancellation of large numbers to yield relatively 


small ones can be avoided. In the denominator of the amplitude in eq. (27) only the 
sum of the expansion coefficients is constrained, di + fi < Large LECs 

matching this condition pose a numerical problem in the course of diagonalising the 
Hamilton matrix for A > A systems. 

2. The computational cost is significantly reduced because of the scalar operator struc¬ 
ture of the interaction. 


Projecting the potential in partial waves, the LECs Ai can be expressed in terms of spec¬ 


troscopic coefficients (table XI in app. |^. Imposing the constraints in 

table IV, i.e., set = 0 = results in an interaction where the leading-order 

part - given in eq. (22) &: (25) - is still iterated, but the NLO part 
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contributes only perturbatively. The latter is now of the form, 


yNLO 

pert 


^ ('e ^4 + A' 2 ai ■ aj) 


Kj 


+ (^3 + ■ ^j) {e ^4 , v^l + ^A^e ^4 ■ aj) r|.) 


E —— A/ -> -> 

e 4 o'e 4 jkA^^^Ti-Tj 


( 29 ) 


i<j<k 

cyclic 


following a prescription given by m- The three coefficients corresponding to momen¬ 
tum independent operators are functions of the two independent NLO LECs and the 
three LO LECs, ^ 3^112 = (^3 4 ; ^ 3 Ni,i, 2 ) • This relation enforces the constraint 
hm T^^{E) = ^im which is not satisfied automatically if a cutoff regulator is 

used. Explicitly, in the limit of taking the energy to zero, the NLO calculation has to 
yield the LO values for the scattering lengths and binding energies of the deuteron and 
triton. This is guaranteed, first, by an expansion of the LECs in the low-momentum scale 
in addition to that of the operators. 


A. = ^ Af' (A) , 


(30) 


k=-l 


where A parameterizes the regularization scheme. Here, this expansion is resembled in the 
A[ 2 terms. Second, one has to match the LECs to an appropriately expanded amplitude, 
where my choice is listed below. 

Prior to this work, the strictly perturbative treatment of EET interactions was sucessfully 
applied to three and four fermions in a harmonic trapES using short-range interactions simi¬ 
lar to the ones here. To results for three-nucleon scattering, analyzed with this approaches], 

Eurthermore, calculations SSI in a discretized space (lattice) were 


we compare in ch. 


VIC 


conducted for ground-state energies of A = 3,4, 6 ,12 nuclei with a chiral effective field 
theory which also treats NLO corrections perturbatively. 

Eor this work, the interaction is derived for the analysis of the 5- and 6 -helium system, 
and proceeds in two steps. The LO LECs are determined to yield the unperturbed states 
which are then used to calculate the NLO matrix elements: 
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TABLE VI: Numerical values of the cutoff A and the LECs of the pNLO EET('jt) potentials. 


A 

[MeV] 

V 

^LO ^LO 

87r2 87r2 

[MeV] [MeV] 

A'l 

[MeV-fm^] [MeV-fm^] [MeV-fm“^] 

A 1 

[MeV-fm“^] 

A 3 N 1 

[MeV] 

300 

as , B{d) 

-81.480 -6.9786 

-43.198 

-60.593 

-10.062 

-7.1631 

-3.0509 

400 

as , B{d) 

-136.66 -9.1401 

-129.41 

-190.07 

-18.536 

-12.634 

2.0703 

500 

as , B{d) 

-81.505 -6.9786 

22.954 

9.6193 

1.2529 

0.83394 

-35.996 


1. Iteration of LO vertices is encoded in two- and three-body eigenstates with 
V = {as,B{d),B{t)} and scattering states with E = 0.4 keV: 

: (f +IV’) = ^ 

'-V-" 

=HhO 

(73)^1 : = B{t)\ijj) —)• triton ground state |t) 

^ = E\^) B{a,LO),\a) (31) 

The ground state wave function and the energy of the four-nucleon system are pre¬ 
dictions at LO. 

2. Leading-order wave functions are perturbed by single NLO interactions with LECs 
A'l 2 3 NI 3 4 fixed by the constraints: 

(i) (4“) (15-0)11/^^014^^ ('So) ) = 0 

(ii) (4“) (3^1) 14 ^^ 044 ^ (^^Si) ) = ^ 

(iii) = o 

(iv) = 0 

(v) (a|l/4(;0|a> = -i?(a,LO) + B(a) 

This means using the np scattering lengths, and the binding energies of the deuteron, 
triton, and o. as input. We will also use interactions which fit to B(a) at LO. For 
those, constraints (iv) and (v) are changed accordingly. 

The reasoning behind the unusual choice of input data to fit the NLO LECs is given in 
app. 1^ after technical details of the numerical method and findings in the n-d scattering 
system, necessary for the understanding, have been discussed. Here, we stress that a four- 
body datum instead of an effective range is used to determine NLO LECs. The second 


j ^ ~ scattering states 44^) 

[ E = —B{d) : deuteron ground state |d) 
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condition inclndes the scattering length in the deuteron channel in the input data. At LO, 
B{d) is used to determine the LEG in this channel, resulting in a scattering length 
of a 4 (LO) = . Through condition (ii), the experimental at is reproduced at 

pNLO. We derived three potentials (table |VI[ ) and applied them to the five- and six-body 
system (fig. & [Tt]). All calculations so far, if not stated otherwise, utilize the variational 
method which I will introduce in the following chapter. 

IV. THE RESONATING GROUP METHOD 

For the approximate solution of the few-body problem, I chose a variational approach 
in which the total wave function is sought as a linear combination of products of bound 
subsystems. This decomposition into fragments is based on the idea of nucleons jumping 
resonantly from one defined group within a nucleus to the otherThis ansatz is 
expected intuitively to be particularly handy for systems where the jumps occur in short 
time intervals relative to the resting time in a certain configuration. If in addition the state 
of interest has one dominating component, in which it spends most of its life, this approach 
allows for an efficient yet accurate description. For this work, the method was chosen for 
its following properties: 

• Calculation of nuclear bound- and scattering states in principle for any number of 
particles A in one framework. 

• Direct access to an approximation of the totally antisymmetric wave function is pro¬ 
vided. 

• Individual basis states can be classified as describing high- or low-energy modes, and 
hence renormalization arguments provide a guideline for a systematic refinement of 
the model spaces. 

• Correct treatment of long-range forces, e.g., Coulomb interactions between charged 
constituents. 

• Multi-channel framework allows the treatment of breakup reactions (limited to two- 
fragment asymptotic states). 
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I will review the calculation of observables in this framework, while details about the eval¬ 
uation of coordinate space matrix elements and model spaces are dealt with in app. 0 and 
during the discussion of the results for the respective systems. 

In our incarnationSi] , the Resonating Group Method (RGM) treats the scattering of 
two - in general composite - nuclei and decomposes the translationally invariant Hamilton 
operator - the center of mass motion is separated off - accordingly into fragment-internal, 
short- and long-range inter-fragment parts (relate to fig. for an example and notations): 

^ + ^f2 + T;el + Vboul,rel + - Vboul,rel • (32) 

'-V-' 

= Rhort 

This motivates a factorization of relative motion from the channel function in the ansatz 
for the total wave function: 

JV(ch) j[= ^ (-l)^R 

(33) 

A^(ch) : number of channels 

The wave function is antisymmetrized (M) with a symmetric group S{A) for A particles 
decomposed into double cosets 0^ . This decomposition reduces the number of spacial matrix 
elements to be calculated to one representative for each double coset instead of one for each 
permutation V. The superscript m labels the boundary condition imposed on the relative 
motion wave function: 


^ ^ ^muQn ^ ^ 

V 

with 


R OC fcm(fl) - rcm(f2) 

^]^rnn . variational parameters 


, (34) 


hm = (^fi + kik iT + , (35) 

\R\^oo ^ > 

i.e, {Gk is the regular (irregular) Goulomb function®) for the relative motion of the 
two fragments with a kinetic energy ^ . In its current implementation, the RGM 

imposes the boundary condition 

ui{kR) = (hYikR) - e2*'5*/i+(feR)) , with hf{kR) = Gi{kR) ± iFi{kR) , 

(36) 

to adopt the convention of [H] for the scattering matrix. The boundary condition allows 
regular waves in channel n only, while admixtures of other channels, and the corrections for 
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mn 
1 / ’ 


relative separation \R\ of the fragments where 14hort(^) / Oj are encoded in a™" and U 
respectively. The variational parameters are determined by the extremal condition (Kohn 
variation®] ): 


5 = 0 


(37) 


In general, the variational approach yields an asymmetric reactance matrix, an artefact 

eliminated by a so called Kato correcture®, after which the S-matrix is derived: 

^ _ 1 + ia 
1 — ia 

The divergence of Gi{p) at /? = 0 is regularized here by using 


(38) 


G^kR) = £ 


yj=2l{n)+l 


I • G^iR) 


(39) 


instead of Gf in a partial wave expanded eq. (34). This approximation, although negligible 
in the determination of a™'”', has a significant effect on matrix elements involved in the fit 
of perturbative NLO LECs. Specifically, limA^oo (V’o ^ ||) receives no contribution 
from G for values of the regulator parameter /3 that suppress the function over the range of 
support of the potential (see also app. 0. Increasing (3 to mitigate the effect of the regulator 
for potentials of about the same range requires in turn a larger model space to approximate 
the G. Therefore, I choose relatively small cutoff values A corresponding to longer-range 
potentials, and used available independent bound-state data to determine the LECs for the 


potential of eq. (|29|) 


Like the Coulomb functions as solutions to the relative motion part, the channel wave 
function - containing the bound states of possibly composite fragments - is expanded in 
a finite basis in order to render the method practical. Our implementation employs a 
Gaufiian basis. Another approach, recently developed®], uses a harmonic-oszillator basis 
and therefore orthogonal basis vectors - an advantage over the GauBian basis. I distinguish 
between three types of expansions: 

1. Cluster expansion: Of all possible groupings G, only a subset Q' is considered in the 
ansatz, e.g., from 

^(®He) = {^He-2n, t-t,^ He-3n, d-d-2n, (pp)-2n-2n} 

the 6-helium ground state was found®] to be dominated by a single configuration 
Q' = {^He-2n} (see fig. 


16 


for a sketch). In a scattering calculation at energies below 
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a certain threshold, closed channels contribute in the form of so called distortion 
channels. Such channels do not have an asymptotic tail and are used to enlarge the 
basis for a better approximation of the wave function in the region of configuration 
space where the colliding fragments interact. Therefore, those channels are included 
to increase to variational space although they do not correspond to, e.g., possible 
breakup of a composed system. 

2 . Partial wave expansion: Orbital angular momenta up to ^m,rei(int) are considered for 
the relative (fragment-internal) wave functions. For the bound states of the fragments, 
I include all coupling schemes of fragment-internal angular momenta I consistent with 
its total orbital angular momentum L. For an ^4-nucleon fragment there are A —1 Ts 
and A — 2 intermediate couplings I' to specify. Likewise, all coupling schemes of the 
individual spins of the fragment’s constituents to yield a good total spin quantum 
number S, that can be coupled with L to good J, are included in the wave function 
ansatz. Only the A —2 intermediate spin coupling values have to be specified because 
of the conserved nucleon spin of ^. 

3. Gaufiianization: The radial parts of the partial-wave-expanded functions are ex¬ 
panded in a Gaufiian basis parameterized by a finite set of width parameters, 7 . 

The regularization parameter ArgMj introduced in sec.jnj is then a function of: Q', a set of 
possibly different /max; and 7 . The convergence of the calculation was assessed with respect 
to all three expansions. The decomposition in two fragments, and the angular momentum 
coupling scheme are manifest in the ansatz for the channel wave function: 


V’ch 




(g) 



where the abreviation for the coupling of two spherical tensors. 


(40) 


[5™^ ® , (41) 

ms,I 

is used. The channel spin Sc is coupled from the total spins Jfi of the two fragments 
whose wave functions are denoted . Appropriate labels for components of the spherical 
tensors are understood. The implied coupling scheme - relative orbital angular momentum 
Lrei with Sc to a total J - differs from the fragment-internal scheme for Jji. The ansatz for 


30 







FIG. 3: Example of a decomposition of berillium-8 in two helium-4 cluster. Subscripts indicate the 
number of the proton or neutron used to label spin- and orbital coordinates. Within a fragment, 
standard Jacobi coordinates are used: Pi = fi — . 

the fragment wave function is: 


E E 


/es' ij.sj 
j 


A-1 

[n- 

n=l 




(Pn 


1 L{j) 


n J 




with 


—1} , n€{l,A — 2} 

h = {{Uj),lUj)),L{j)} 


and 7 , = { 7 ™(j)} 
Vl (f) = |r|W; 


( 42 ) 


The two fragments of the example in fig. [^also represent two different values of /, i.e., 
substructures within a fragment. The superposition coefficients c are determined by solving 
the generalized eigenvalue problem because of the nonorthonormal nature of the basis: 


{ML I H I A(f)i ) = E I Acl)i) , with m,n = {/, j, Ij, Sj, 7 j} defining a basis vector. 

'-V-" 

in gen. 

=N ^ 1 

( 43 ) 

The eigensystem is solved by diagonalizing H' defined as 


N-'2HN-{ -E1 = 0 . ( 44 ) 

H' 

The inversion of the norm matrix N increases the numerical uncertainty and renders the 
method useless if N has eigenvalues of order 10“^^ or smaller. Eigenvalues of this magnitude 
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can result from basis states which are linearly depended up to differences which cannot be 
resolved by the numerical precision, and have to be excluded. In app. 1 a three-nucleon 
system is considered as an example for the construction of the fragment wave function 


eq. (42) and the following determination of the ground-state wave function. 


With fi and gi expanded in Gaufiians - in practice, a basis parameterized by one of 


the sets 'jj for the states in eq. (42) is used - the matrix elements necessary to determine 
the reactance matrix a™'"' in the channel basis require a basis transformation resembling 
the different angular momentum coupling schemes and hence can be written for a generic 
operator 6 as a linear combination of elements of the form: 


Mmn = {4’m^ (notation def. in eq. (|4^). 




(45) 


The spacetime symmetry of the theory (ch. B allows only spherical tensor operators o(r, rrir ) 
of rank r = 0 which makes the following separation in coordinate- (do) and spin part (d<j) 
possible: 


H0>0) = '^Oo{q, 


rrin 


ds{q,mq) 


JzllL 

y/^q -\- 1 


(- 1 )^ 


(46) 


={qmqq-mq\00) 

The operator structure for the EFT potentials in eq. ( [2^ &: ( |25[ ) matches this form if the 
tensor operator is transformed, Ui ■ raj • r — ■ aj = \/5[[(Ti (8) aj]'^ (8 [r (8) r]^]°. Using 

Racah algebraSS, the matrix element can be written in the form 


Mmn = Sjj'SmM' ^ (“ 1)^C* (j)c(j) (-1)^ 
VeSiA) 


L+2S+S'-J I S L j\ , 


L'S' 


{L'\\Vdo{r)\\L){S'\\Vd,{r)\\S) . 


The Wigner-Eckart theorem defines the reduced matrix elements 

(5'||Rd.(r)||5) = {S'S'\Vds{T,S' - S)\SS) 


(47) 


(48) 


The magnetic quantum labels are conventionally chosen to be maximal. Eq. 48 holds 
analogously for the reduced coordinate space matrix element. 
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Evaluation of the spatial matrix elements. 

coordinate-space matrix elements necessary to obtain M^n 
(triplet) of interacting particles ij{k) is: 


The basic structure of all 


in eq. (45) for a given pair 


/ dsi...dsA-ie 




' ■Wyir.ml^e^^Xvy, 

\ At=l 


n=l 


= Qn 


(49) 


This expression is valid if the following conditions are satisfied: 


• The single particle coordinates ri are transformed in Jacobi coordinates Si via an 
orthogonal matrix in bra and ket. 

• The permutation V is applied to the right, he., to the coordinates of the potential 
operators and the ket. 


The part of an operator acting in coordinate space, do, is expressed in terms of 


Gaufiians and spherical harmonics, as in eq. (26). The choice for the regulator re¬ 
sults in Gaufiian radial potential functions, and therefore, no additional expansion is 
necessary. 


The Jacobi coordinates of the bra are to be chosen as the independent variables, i.e., 
all other coordinates are expressed in terms of them. 


These operations define the transformation matrices p and ^ (eq. (49)), both depending 
on the coordinates of the interacting particles and the permutation. The Qn are hence, in 
general, linear combinations of the independent variables Sn- z equals the total number of 
inter-cluster coordinates in bra and ket plus the number of spherical harmonics from the 
potential operator. For example, the operator oc Ag in eq. ( [2^ adds 2 to z. The two major 
simplifications which lead to an analytic expression for J{V) are: 


• the diagonalization of using a transformation matrix Txx' with Txx = 1 and 
T\X' = 0 for A > A', 


expressing the spherical harmonics in terms of their generating functionSS, 

L 

(b-f)^ = b^ ^ CLmb~'^yLm{r) , 
m=—L 


(50) 


with b={l-b\i{l + b^),-2b) , and Clm = {-y) L'-J (2L+l){L-m)\(L+m)\ 
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An example of this prescription is given in app. 0 for the calculation of the spacial part of 
the operator oc Ag - as the only structure not present in previously considered potentials, 


like AV180, CD-BonnI3 - in eq. (26). 


Evaluation of the spin matrix elements. 


r in eq. (26) acting on particles i and j is given as: 


A generic spin operator with rank 


Os{ij-,r) = d^{i) ^ {ls{i),di} ■ (51) 

m 

The constants Cm represent the coupling of spherical tensors of rank 1 {a) and scalars 
(1) to a spherical tensor of rank r. The transition from the spherical components to 
raising and lowering operators is done via with known action on spin 

eigenfunctions. Accordingly, the state \SS) is expanded in a linear combination of products 
of single particle spin eigenfunctions 

|55) = • ... • I^tun) . (52) 

k 

The Ck contain the Clebsch-Gordan coefficients of all couplings. The matrix elements 
{S'S'\'Pds{ij-,r)\SS) are calculated by applying ds{ij,r) to the right, yielding another prod¬ 
uct function. Now all permutations V are determined which result in non-zero matrix 
elements whose value depends on the specific operator coefficients. 


V. A=2 


The two-nucleon system is relevant for this work because, first, it provides most of the 
data used to determine the LECs, and second, the sensitivity of observables to Arqm can 
be assessed by a comparison to predictions by other numerical and analytical calculations. 


The parameter Arqm is defined in ch. IV as a set of fragmentations Q', a set of maximal 
orbital angular momenta with corresponding coupling schemes, and a set of width param¬ 
eters covering an interval which is cut off at short and long distances. Although Arqm is a 
renormalization parameter in the same sense as the regulator width A, probing for incon¬ 
sistencies in the power counting is less practical with it because it is not straightforward 
how to change Argm such that short-distance physics only is affected, especially for A > 3 
systems. The theoretical uncertainty in an observable at order n must not significantly 
exceed , i.e., VA > M predictions have to agree within this error margin. For Arqm, 
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there are corresponding inequalities, one for the angular momenta Im, and another for the 
width parameters (assuming G' = G). Analogous to the A dependence, the soundness of the 
power counting guarantees that under a modification of the model space, i.e., an addition or 


removal of basis states of type eq. (42) with ln{j) > ^max and 7 n(i) ^ 7 max - corresponding 

min 

to high energy modes - predictions for low-energy observables change only at order (^ 


Choosing a complete basis is equivalent to A —?■ oo. For two reasons I favor to operate in a 
numerically complete basis, which does not exclude short-distance modes from the start: 


In practice it is more convenient to vary A once a model space with Arqm is chosen. 


• It is unclear how a certain ArgMj imposed in the A = 2 sector, is consistently trans¬ 
lated to A > 2 systems, e.g., what triton model space is to be used given a deuteron 
model space to yield the same Arqm? Therefore, a change in the renormalization 
scheme is expected from A = 2 to any larger system. Harmless in the case of an 
approximately complete basis, a smaller model space could result in a significant de¬ 
pendence of the LECs on Argm in addition to that on A. Hence, LECs determined 
from A = 2, 3 data and used to predict A = 4 observables would falsely suggest a 
failure of the power counting because of the change in Argm- 


The two-nucleon model space, as specified in table |VH[ consists of a set of 20 width pa¬ 
rameters and a pure S-wave (/i = 0) for the LO and perturbative NLO calculation, while 
admixtures of D-waves {li = 2) are possible in the nonperturbative NLO approach. The 


small values for 71 are included to expand asymptotic states hf (eq. (36)) while the large 
widths are needed for the bound state wave functions, for which structure at distances oc A“^ 
is relevant. The LO wave function is subject to two boundary conditions: a rapid decrease 
towards zero, r —)• 0 , starting at |r| ~ A“^ and a certain asymptotic decay for |r| ^ M~^, 
independent of A. With a basis like tcl20, the second condition can be met for binding 
energies of the order of the physical B{d), but any such basis, with a maximal width 7max, 
will not be able to simultaneously reproduce the falloff for A > ^ 7 max (eq. ([M|). This last 
assertion is verified by checking the consistent prediction for the deuteron binding energy 
in two calculations: 
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TABLE VII: Two-nucleon RGM basis used for the determination of LECs and predictions. The 
values were optimized for a potential model calculation and the basis is labeled wl20. wl2 refers to 


a basis with widths 71 (j) = 0.1 • 71 (j)- The notation is defined in eq. (42). The deuteron state is 
specified by spin s = 1, a singlet deuteron by s = 0, while the neutron-neutron and proton-proton 

I 

systems have s = 0, too. The total isospin is fixed by: I + s + t = odd. 



widths [fm 



dim 

LO,pNLO 

7 i 0 ') £ W120 = < 

f ' 

129.567, 51.3467, 29.4729,13.42339,8.214456, 

4.447413, 2.939,1.6901745,1.185236, 0.84300, 

:> 

0 

II 

20 

nNLO 

71 O’) ^ W120 

0.50011,0.257369, 0.13852, 0.071429, 0.03852, 

0.01857,0.009726, 0.005619, 0.00277, 0.00101 
>. / 

r 0 ; i < 20 
h{j) = 1 

^ 2 ; i > 20 

40 


1. The LO LEG in the deuteron channel, Ct = C\^+€ 2 ^, is htted to at with a standard 
numerical integration, using Mathematica and not the RGM, of the Schrodinger 
equation for M < A —?■ 00 . 

2. With this Ct and the RGM model spaces wl2 and wl20, an RGM deuteron binding 
energy B{d,A) is calculated. 


Figure 1^ shows consistent predictions for B{d) in wl2 and tcl20 up to a threshold value 
1.4 GeV. For A > u;12 fails in expanding a bound state, while rcl20 does so 

above ~ 3.6 GeV. Furthermore, the predicted B{d) converges as expected for A —>■ 00 


to the analytical prediction given by eq. (23). 

The accurate expansion of scattering states for leading-order potentials and the con¬ 
verged variational determination of scattering phase shifts is demonstrated by comparison 
to predictions by a direct integration of the Schrodinger equation with Mathematica for the 
neutron-proton ^S'l-phase shifts in hg. 

I conclude from these analyes that the dependence of two-body LEGs, obtained with a 
basis t(;120 for A < 2.5 GeV, on Arqm is insignihcant relative to the systematic uncertainty 


of a NLO calculation with EFT(^). With this basis and the limitation on the cutoff, the 
RGM is an accurate numerical method for A = 2. For the next larger system, the triton. 
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FIG. 4: Binding energy B{d) of the deuteron in leading-order EFT('^) as a function of the Gaufiian 
momentum cutoff A for two different RGM model spaces. The space wl20 contains narrower width 
parameters, and hence can be used for the potentials corresponding to larger values of A. Using 
the triplet neutron-proton scattering length at as experimental input, the effective range formula 
predicts B{d) = 1.41 MeV. 



FIG. 5: Relative deviation between the np phase shift results obtained with the RGM and a 
numerical integration of the Schrodinger equation. A leading-order potential with one parameter fit 
to at and a cutoff A = 1.5 GeV was used. Notice the scale on the ordinate. 

the basis wl20 is used to expand the deuteron fragment and a part of the width set W 120 
to expand the radial function of the second neutron’s coordinate in order to translate this 
accuracy to the A = 3 system. 
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TABLE VIII: Triton and S-helium RGM model spaces. Quantum numbers, which are not specified, 

run over all values compatible with the total angular momentum and parity, and coupling 

rules, e.g., the li ^2 = 2 set contains basis vectors with total orbital angular momentum L = 0,1,2. 

fl 2) 

The nNLO basis includes all possible combinations of /i, ^ 2 , and s. The sets W 120 subsets of 
W 120 containing every other value starting with the first (second) width. 



h 

I 2 

S 

width sets 

LO.pNLO 

0 

0 

1 

72 (i) e >Vi2o,7i(i) e 


0 

0 

0 

72 (i) € Wi2o,7i(i) e 

nNLO 

{0,1,2} {0,1,2} {0,1} 

71,2 given in [33] 


dim 


200 

200 

220 


VI. A=3 


While two-nucleon observables at momenta p <C are parameterized accurately by 
the two scattering lengths only, an additional parameter is needed for the three-nucleon 
amplitude (see discussion of eq. (EH))- The dependence of observables on the renormalization 


parameter - manifest, e.g., in the Phillips line, the universal ratio eq. (19), or the Thomas 
effect - expresses this feature of the nuclear interaction. The estimate for the expansion 


parameter in eq. (16) would be false if a similar promotion of a higher dimensional operator 
occurs - i.e., a LEG violating eq. ([^ with <5 > 8. The cutoff dependence of three three-body 
observables is investigated here at LO and NLO. The single additional parameter allows 
for a tangible graphical representation of this dependence by plotting the prediction for an 
observable as a function of another for varying renormalization parameters. A renormalized 
order n interaction corresponds to a locus of points bounded by a rectangle around the 
datum. The size of this bounding box resembles the theoretical uncertainty with respect 
to different unresolved microscopic structures and relating these dimensions at successive 
orders assesses the expansion parameter of the effective field theory. This is how the naive 


estimate of eq. (16) is tested here. All three-body observables were obtained with the RGM 
with variational basis defined in table VIII In addition, we show in app.j^with the example 
of a smaller basis explicit formulas for the triton trial wave function. 
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2.8 



B(t) [MeV] 


FIG. 6: The correlation between the triton charge radius and its binding energy. The narrower 
band mapped out by the EFT(/) NLO results from the RGM (blue squares, black dashed lines) 
as described in the text, compared to our RGM LO results (gray shaded area) and a LO calcula¬ 
tion by [5D] (dashed gray lines), demonstrates convergence, also to the datum (exp. value) and a 
customary finite-nucleon-size-correction of it (pp). The values from AV18(+UIX) (black triangles, 
RGM calculation for this work) and a variety of other potential models (gray diamonds [51] |S2] ) are 
consistent with our proposed EFT(]t) NLO correlation band. 

A. Triton charge radius 

The variational RGM approximates the ground-state wave function {f\t) of the triton as 
the eigenfunction corresponding to B(t). Hence, a numerically complete variational basis 
and therefore a converged B{t) transients into convergence for all other observables derived 
from the ground-state wave function. In particular, we consider the average charge radius 
which is calculated from the wave function using: 

Yj 

where f) is the position and (1 -|-r?) the charge operator of the z-th nucleon. An analysis of 
the correlation between B{t) and {r‘^)Y, as motivated above, will provide information about 
the convergence properties of EFT(^), and is presented in fig. Predictions are shown for 
the LO and nNLO EFT(^) interactions of this work, a LO EFT(^) calculationextracting 
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observables by solving Faddeev integral equations, various nuclear force modelsEH^, and 
two RGM calculations for the AVIS® two-nucleon- and the AV18+UIX® two- and three- 
nucleon interaction model. For the discussion of this graph, I define two terms: 

correlation line: {{B{t) (Csni) , o (Gsni) ) , —oo < Csni < oo } , 

i.e., fixed A, while Arqm needs to be adjusted for shallow or deep bound states 
relative to the state at Csni = 0. In this section, is considered as the correlated 

observable o. 

order n correlation band: {( B{t) (A, Arqm, C'sni, B') , o (A, Arqm, Qn!, B'))] , 

i.e., all renormalization parameters A,ArgmjC' 3 NI and conditions at n-th order, V, 
accessible within this RGM approach are varied. 

The EFT formalism identifies the band as the area in which the predictions of any two- 

body interaction - parameterized by as^t or as,B{d) - are predicted to be. In that sense, 

the bands are the universal three-body consequences of specific two-body interactions. By 

definition, the bands are unbounded, resembling the nonperturbative dependence on the 

renormalization parameters. Renormalizing the theory with an additional constraint on the 

three-body spectrum, enforcing a “triton” energy in an interval around the datum of width 

set by the uncertainty of the considered order, automatically constrains all other low-energy 

observables derived from the wave function corresponding to this bound state, or scattering 

observables with one nucleon impinging on the deuteron. 

Turning now to the results in fig. one first confirms the expected behavior of an 

increasing for more loosely bound systems. Results of two EFT(y') LO calculations 

are shown. The RGM LO potentials map out a band (gray shaded area) which includes the 

1 

datum. In addition to the experimental value, — (Rp) — 

is given as a measure which takes into account the finite extension of the nucleonsES, 
(Rp) = 0.769(12) fm^ and (R^) = —0.1161(22) fm^, plus the Darwin-Foldy termES (red pp 
cross in fig. [^. The corrections constitute a partial removal of higher-order effects. Hence, 
it has to be an element of the LO band. The important reference point here is the datum 
itself, because all finite-size and relativistic corrections are expected to be included order 
by order in the EFT expansion. 

The lower and upper LO line from [50j (dashed gray lines in fig. result from a fit of 
the LO LEGs C*]"® to either ag^t or B{d), a^, respectively, providing some measure of higher- 
order effects. For the RGM LO calculation, we used ag^t (red solid lines) and B{d),ag (red 
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dotted lines, with the closest solid red line corresponding to the same A) to determine the 
LECs, and varied the cutoff from 200 MeV (top edge) to 1.6 GeV (lower edge) in steps 
of 100 MeV. The dependence on V, i.e., using B{d) or at as input, measured in terms 
of relative shift of the correlation line, becomes comparable to that observed by doubling 
A = 800 MeV to 1.6 GeV but is found negligible for a doubling of A = 200 MeV. We choose 
the lower bound for A to be in the region of the pion mass. The upper bound is set by 
the observation that the band width is increased little by increasing A from 800 MeV to 
1.6 GeV compared to the broadening from 200 MeV to 400 MeV. The band appears to be 
saturated, with the bulk of its width coming from the region A G [200; 800] MeV. For fixed 
A, was varied to generate the correlation lines. 

The quasi-exact Faddeev calculation of m and the RGM results do not overlap. With 
fig- i I demonstrated that numerical inaccuracies of the RGM are negligible, and 
conjecture that the non-overlap can be traced to the differences between the regular¬ 
ization schemes employed. The Faddeev calculation uses the separable cutoff function 
/k(p,p') = exp (—exp (—with k > 1600 MeV, while the RGM uses the non- 
separable regulator of eq. pO] ) plus an implicit regulator imposed by the finite number of 
width parameters. Therefore, both methods combined can be viewed as providing a check 
of residual regularization-scheme dependence. A conservative estimate of LO effects is thus 
the range of results covared by the combination of both methods. The LO accnracy at the 
physical triton binding energy is thus ±0.6 fm, and the measured charge radius happens 
to lie right in the middle of the LO band. Additional regularization schemes, e.g., with a 
range of cutoff values in the Faddeev approach, are thus conjectured to lead to an overlap 
between the two bands. 

At NLO, the potentials without 3NI map out a correlation band which is more narrow 
than its LO counterpart, and again contains the datum. The width of the band can be 
estimated by considering the difference between the results using the extreme cases provided 
by potentials and They produce the same NN scattering lengths and nearly the 
same triton binding energy, and their deuteron binding energies differ by less than 10%, 
while their triton charge radii differ by about 10%, consistent with the expectation of 
a NLO calculation. Finally, they are based on two significantly different cutoff values, 
A(^V^) « 400 MeV and A(®V)f) 650 MeV. We therefore can base an estimate of the NLO 

band on the range mapped out by varying the 3NI for these two potentials. These lines are 
included in fig. 
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Variation of the cutoff A, of the fitting-input P', and of the 3NI lead therefore all to 
similar assessments of the uncertainty of the theory at NLO. At fixed B{t), the charge radius 
varies from LO to NLO by < 30%, in agreement with the power counting which estimates 


the correction according to eq. (16) to be approximately i. 


^ and the observed 


All three values, namely, ±0.2 fm from the naive estimate Q 
convergence from LO to NLO, and ±0.1 fm from the above mentioned difference between 
and would be equally valid estimates for the theoretical uncertainty of this NLO 
calculation. Using the more conservative error estimate and the assumption for the NLO 
error band center at 1.6 fm at the experimental B{t), EFT(^) predicts a value 


(r2)^Y^(LNLO) = (1.6 ±0.2) fm 


(54) 


within the band of the leading-order value LO) = (2.1 ± 0.6) fm as 

quoted from |50j . The NLO value is found in good agreement with experiment®), 
(r^)^/^(L exp) = (1.63 ±0.03) fm. 

Another argument in favor of our definition of the correlation band is provided by the 
results of the two phenomenological models, AVlS^l and AV18+UIX0. In general, EFT 
predicts that the results of a potential which reproduces or shares input observables, at least 
to the accuracy required at the considered order, deviate from the results of an appropriate 
EFT potential by less than the theoretical uncertainty of the EFT values at this order in 
the applicability range of the EFT. This criterion is easily met by AV18(+UIX), and hence 
its predictions have to be consistent with the proposed correlation band. The two-body 
potential AV18, reproducing the Nijmegen phases much more accurately than required to 
fall into this category of potentials, is expected to yield a value at a position within the 
((^^)ch(i) — il(t))-band. The prediction for the triton charge radius of AV18+UIX is expected 
to deviate less than 10% from the experimental datum within the error band, because 
this model has a three-body interaction added to reproduce the experimental B{t). Both 
expectations are consistent with the results shown in fig. 

In conclusion, the results for these three-nucleon observables show that although the 
potentials are approximately NN phase-shift equivalent, they differ in their predictions 
of three-body observables. As mentioned above, a proper renormalization of the theory 
requires therefore one three-nucleon contact interaction. Setting this 3NI to zero in the 
potentials the expected dependency of observables in A > 2 systems on how the 

unobservable short-distance physics is modeled is observed. Different short-distance physics 
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FIG. 7: Binding energy difference between ^He and for EFT(|f) potentials (filled squares) com¬ 
pared to various isospin-invariant potential models calculations!^ (circles) and RGM values for 
AV18(+UIX) (triangles) which contain charge symmetry breaking terms. For the potentials 
(solid) and (dashed), a smooth variation of the 3NI leads to the two correlation lines. The 
upright cross is the experimental value without the contribution from CSB terms Efil. 

is modeled by the potentials not only by varying cutoff values but also by differing sets of 
LECs for the same cutoff, while the scheme- and regulator-dependent three-body interaction 
parameter is chosen as zero. With this parameter fitted to the triton binding energy, the 
prediction for the triton charge radius is consistent with experiment within the expected 
uncertainty range. A significant convergence from LO to NLO is observed. 


B. Trinucleon binding energy splitting 

In fig. the RGM results for the splitting between the binding energies of 3-helium and 
the triton are shown for the potentials At NLO in the pion-less EFT, the strong inter¬ 

action is isoscalar and hence does not break charge symmetry. Therefore, charge symmetry 
breaking (CSB) comes in our RGM calculation only from including the Coulomb interaction 
between the protons in ^He. The EFT(^) results, hence, only show the model-independent 
contribution of Coulomb interactions to the trinucleon binding energy splitting. The cor¬ 
relation band is mapped out by the EFT(^) potentials with zero 3NI and by smoothly 
varying this three-body parameter for potentials Both approaches result in correla- 
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tions consistent with each other and with the results of the other NLO potentials. At the 
experimental triton binding energy, this leads to predicting 

{B{t) - Bi^Re)) = (0.66 ± 0.03) MeV . (55) 

Here, we estimate the theoretical uncertainty by the spread of the phase-equivalent NN 
potentials as in the previous section, e.g. comparing The a priori error estimate 

at NLO of < 10% gives a larger uncertainty of ±0.07 MeV. However, one should keep 
in mind that including isoscalar strong interactions by higher-order terms of the effective- 
range expansion has identical effects on the strong interactions inside the triton and ^He. 
These effects cancel out in the difference and only survive indirectly, as the strength of the 
Coulomb interaction in a system is also correlated to its size. We therefore quote the width 
of the correlation band as the error estimate of our calculation. 

This value deviates by about 0.1 MeV from the experimental value^S of 0.764 MeV. In 
line with the argumentation above, we attribute this difference to isospin-breaking charge- 
independence breaking (CIB) and CSB interactions coming from the explicitly broken chiral 
symmetry in the strong and electro-weak sector from which only the parts resulting in the 
Coulomb force have been considered in this calculation. They enter in EFT(^) only at 
higher order. To support this assertion, results of the potential models AV18(+UIX), which 
contain CSB interactions, are included in fig. Both potentials are not elements of the 
correlation band suggested by the EFT(^) points but agree with a shifted band, centered 
around the datum. In contrast, the values from charge-symmetric potential models in fig.[^ 
lie within the NLO EFT(y') band. This leads us to predict a model-independent CSB/CIB 
contribution to the binding difference in NLO EFT(^) at the experimental triton binding 
energy of 

{B{t) - H(3He))^®®^^^^ = (0.10 T 0.03) MeV , (56) 

anti-correlated with the Coulomb contribution to give the experimentally established dif¬ 
ference. This is to be compared with the contributions from two- and three-nucleon CSB 
interactions which stem from Chiral Effective Field Theory, Breit and vacuum polarization 
corrections, and from corrections to the kinetic energy operator. These were calculated^SSl 
to sum up to (0.112 =F 0.022) MeV, leaving about (0.652 ± 0.022) MeV for the soft pho¬ 
ton effects, dominated by the Coulomb interaction. This is in perfect agreement with the 
EFT (if) result. 
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C. Neutron-deuteron scattering 


The following analysis of the neutron-deuteron (n-d) scattering system in the triton (total 
spin-^) channel serves as the first example of a scattering calculation in the RGM framework 
with leading-order EFT(y') potentials. The RGM results and their comparison to previous 
calculations motivate the peculiar choice for the input data V for the determination of the 


pNLO coefficients (see ch. III). 

The effective range expansion, as introduced in eq. expands the amplitude around 
zero momentum and yields a deuteron binding energy which differs by about 0.6 MeV 
from the datum at LO. Effective range corrections to B{d) are then large relative to the 
ones in an expansion around the deuteron binding momentum 74. The latter expansion, by 
construction, reproduces the experimental B{d) already at LO. In [58] it was shown how the 
convergence rate of EFT(^) renormalized by a matching to the Z-parameterized two-nucleon 
amplitude is improved for three-nucleon scattering observables compared to an equally 
admissible theory with LECs fitted to the ERE amplitude. The former parameterization 
uses the deuteron pole position 7^ and its residue as input, while the latter employs at and 
rt- The intuitive expectation that accuracy in the reproduction of the tail of the deuteron 
wave function transients to a better description of low-energy elastic neutron scattering off 
the deuteron was confirmed^^ and is demonstrated by the faster convergence of EFT(^) 
from LO to NLO. Naively, I expect a faster convergence in the ®He system with all of its 
subsystems bound close to data in analogy. Following this intuition, B{d,t,a) are used as 
input and not, e.g., effective ranges from the two-nucleon sector. The improved accuracy 
is also expected in low-energy a-n scattering if R(a,exp), i.e., the target, is reproduced 
accurately at LO with respect to its binding energy. Important differences between the n-d 
and n-a reactions which could account for a different result are: 


• The a state is bound deeply below the lowest breakup threshold compared to the 
deuteron, and relative to the low-momentum scale of EFT(^). Neutron scattering on 
an a target should therefore be independent of breakup effects over a larger range of 
energies. 

• The a’s coordinate-space wave function is not totally symmetric due to a reduction 
of the RGM basis, while the model space for the n-d reaction was sufficiently large 
to assume the symmetry of the deuteron’s wave function. 
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In this section, the findings of |39l |58] for the three-nucleon system are confirmed in the 
sense that LO EFT (it) predictions for n-d phase shifts with the RGM are found to be more 
accurate, i.e., less sensitive to cutoff and input variations, for interactions that reproduce 
B (d, exp). Specifically, this is shown for the correlation between the binding energy and 
the n-d scattering length in the triton channel - this is known as the Phillips line and was 
discovered by [59] and explained by [60| - and the n-d S-wave phase shifts below the deuteron 
breakup threshold. The analysis extends the existing EFT(it) (N)LO calculations insofar 
as it provides a more detailed assessment of the uncertainty at LO. In contrast to [58j . the 
LECs are determined here from the ERE as given in eq. (§, ^.e., an expansion around zero 
momentum. 

First, I discuss the n-d phase shifts in the spin-^ channel shown in fig. The 
RGM model space uses the u)120 basis (VII) for the deuteron fragment and the singlet- 
deuteron and nn distortion channels. The width set W 120 defines the basis for the radial 
fragment-relative function. RGM results are displayed for LO EFT(if) potentials with 
200 MeV < A < 1.6 GeV (blue and red lines). Two sets of input data V are used 


for each A to determine the three LO LEGs (eq. (22) &: (25)): = {as^t, B{t)} and 

1^2 = {us, B{d), B{t)}. Values for B{d) of potentials fitted to range from ~ 2.5 MeV 
(A = 200 MeV) to ~ 1.5 MeV (A = 1600 MeV). The predicted phase shifts as shown in 
fig-IH are found much more accurate for potentials whose LEGs are either determined via 
P 2 or their A is such that the experimental B{d) is approximately reproduced even if 
serves as input data. 

In general, the uncertainty at a given order is assessed by the sensitivity of predictions to 
modification of unobservable short-distance physics. Here, if a change in the short-distance 
structure is induced solely by a change in A, the uncertainty in (5„d, measured by the width 
of the band defined by phases whose corresponding potentials use B{d) as input (red dashed 
lines), is smaller compared to the uncertainty found if short distance structure is also varied 
by using differrent input data, namely T>[. The phases (5nd deviate most from the phase- 
shift analysis data (red squares) for potentials which use and whose cutoff is large, i.e., 
A = 800,1600 MeV (blue lines). The LO uncertainty of EFT(it) is estimated by the width 
of the band bounded by the predictions of the potential with A = 1600 MeV (lower blue 
line) and that with A = 200 MeV (top solid red line). This band is not expected to widen 
any further, because, first, the maximal deviation from B{d,exp) is almost reached with 
A = 1.6 GeV (see discussion of fig. [^. Second, lowering A below 200 MeV will cut off not 
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FIG. 8: Results for the Doublet S-wave phase shift (5nd for neutron-deuteron scattering expressed in 
k cot^ (5nd with center of mass momentum k. The vertical line marks the deuteron breakup threshold 
at the LO B{d). The RGM LO predictions (red and blue lines) are compared to a phase shift 
analysisl^ and LO (NLO) momentum space calculationsES (gray (black) line). Predictions shown 
as solid lines use at instead of B{d) (dashed) as input. The blue lines highlight the specific RGM 
results for A = 800,1600 MeV. 

only unobservable high-energy modes but increase the distortion of observable low-energy 
modes. The instability of the phases predicted by the RGM in hg.|^for energies > 1.5 MeV 
resemble the crossing of the deuteron breakup threshold. The position of this threshold 
is set by B{d) and therefore depends on A and Ph For the potential with A = 1.6 GeV 
and Vi, for instance, the collision becomes inelastic at Ecm ^ 1-5 MeV (vertical line in 
fig- [§ compared to a breakup of the deuteron for Ecm ^ 2.6 MeV for LO interactions 
with A = 200 MeV and in EFT(]t). Threshold positions are variable within uncertainty 
margins of a given order. The values of A and the input set V, for example, constrain 
the deuteron binding energy to 1.4 MeV < B{d) < 2.7 MeV. This fact can be utilized as 
follows: 
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The numerical effort for multichannel scattering calculations can be reduced by widening 
the threshold separation without simultaneously losing accuracy in phase shifts at energies 
below the lowest breakup threshold. I consider the four-body problem of the ^He O'*'- 
channel as an example. In this system, the two-channel calculation, with only t-p and 
^He-n open, is expected to be predicted accurately like the n-d reaction of this chapter if 
the renormalization scheme is chosen wisely. This choice is conjectured to set the scheme 
parameters A and V, such that, on the one hand, the two lowest thresholds - t-p and ^He-n 
- are accurately reproduced, e. 5 ., by including B{t) in V, and on the other hand, the d-d 
threshold is lowered, e.g., using A = 1.6 GeV and Pj. The benefit of the larger gap between 
the ^He-n and d-d thresholds is a reduction of the basis size. 

Concluding the discussion of the RGM uncertainty band calculated at LO in EFT(^) 
contains the (N)LO predictions of [39] ((black) gray solid in fig. as well as the partial 
wave analysis data^l. 

By considering the n-d scattering length in the triton channel correlated to B{t), I 
confirm the conclusions drawn from the analysis of Furthermore, the RGM uncertainty 
assessment derived from a change in the cutoff A and input data, can be compared to |58j 
at LO in FFT(^), where only the regulator parameter was varied. The results are displayed 
in fig. where was derived from the phase shift at Ecra = 4 keV, calculated with the 
RGM and the same LO FFT(^) potentials used for the d^d predictions. The correlation 
lines were obtained by fitting the two-body LFCs to ^( 1)2 \ ((blue) red lines) with 

fixed A and a variation of the three-body parameter A 3 N 1 . 

As for (5„d) the uncertainty assessed by potentials based on V '2 with their precise reproduc¬ 
tion of B{d) is smaller compared to that which includes predictions. In more detail, the 
width of the band mapped out by 'D '2 interactions alone (red lines) resembles the sensitivity 
to short distance physics parameterized by a regulator value of 200 MeV < A < 1.6 GeV. 
Hence, a sole change in A would suggest a LO FFT(^) uncertainty of 

^a^d = 0.67 ±0.19 fm , (57) 

with A 3 N 1 fitted to B{t). A further probe into the short-distance structure of the the¬ 
ory, by including interactions using at instead of B{d) as input iV'i), maps out a broader 
band (hatched area in fig. with its (upper) lower bound set by a potential with 
A = (200)1600 MeV and LFCs fitted to V'^ \ {B{t)}), suggesting the less accurate LO 
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FIG. 9: Correlation between B{t) and the doublet S-wave scattering length ^Ond for the neutron- 
deuteron collision (Phillips line). Experimental values are shown as solid black lines. The (gray) 
blue shaded area is the (N)LO EFT(]t) uncertainty band from [^. The hatched area is the LO 
band obtained here with the RGM. (Blue) Red correlation lines mark RGM LO EFT(/) predictions 
for cutoffs 200 MeV < A < 1.6 GeV and input. 

prediction: 

2a„d = 0.67 ± 1.5 fm . (58) 

Increasing A beyond 1.6 GeV is not expected to widen the band significantly, because of 
the relatively small width increase induced by increasing A from 800 MeV (lower solid 
blue line in fig. to 1600 MeV (lower edge of hatched area), compared to the increase 
from 400 MeV (upper blue solid line) to 800 MeV. We have thus demonstrated how the 
limitation to a certain set of input data reduces the uncertainty of an EFT significantly. 
A cutoff variation as the only probe to dependence on high-energy modes does not reveal 
the overall uncertainty which is found considerably larger, i.e., the width of the hatched 
area, representing the correlation for A and D', is wider than the band mapped out by the 
potentials using (red lines). 

For the range of triton binding energies considered here, the RGM LO EFT(y') correlation 
band overlaps with the (N)LO bands ((gray) blue shaded in fig. calculated also with 
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EFT(^) but using a Z-parameterized amplitude to fit LECs and which were obtained by 
calculating the three-body observables from Faddeev integral equations solved in momentum 
spaced. The broader RGM band is explained by the more detailed probe of short-distance 
structure. While in [58] the band spread is induced by a cutoff variation only, for the RGM 
band, input data is varied additionally. Taking the overlap of both calculations as the LO 
correlation band, all NLO predictions are elements of this broader LO band - resolving this 
issue arising in [58] . 


D. Conclusions from A=3 

The analysis of correlated observables is a useful tool worthy to be applied in larger 
systems because: 

• The shape of the area mapped out by varying the renormalization scheme assesses a 
bound for the expansion parameter(s) of the theory. An unbounded area indicates 
an inconsistency in the power counting. 

• A discrimination between insufficient interaction models and numerical effects can 
be made. Gonsider, for instance, a potential tuned accurately to two-body data but 
missing the triton binding energy significantly {e.g., AV18 fig. [^. Also, the position 
and width of a correlation band for a scattering observable - not one derived from 
a bound state, as then, a different, in most cases more tedious, numerical algorithm 
is required - is known from a relatively simple EFT(^) LO calculation. A prediction 
of the model inconsistent with this band is then only related to numerical uncertain¬ 
ties. This understanding will be employed in the four-nucleon scattering calculations 
presented below. 

• A guideline for theoretical model building is given. A modification of an interaction 
to improve the prediction for a selected low-energy observable must simultaneously 
impact correlated observables in a way consistent with the accuracy of the interde¬ 
pendence. 

• Gonflicting measurements - nonoverlapping error intervals - can be categorized by 
falling in or out of the correlation band at sufficiently high order. A comparison of 
the limitations of the theory, e.g., neglected weak interactions or three-body break-up 
channels, and of the experimental setup, suggest an interpretation of the discrepancy 
in the data. 
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The RGM results presented here confirmed for the n-d scattering length and the phase 
shift in the S-wave spin-| channel how a certain choice of NN data for the determination 
of the LECs reduces LO uncertainty. The precisely reproduced deuteron is to be favored 
over, e.g., an accurate fit of neutron-proton scattering phase shifts at Ecm —^ 0 MeV, if the 
LO uncertainty for n-d scattering observables is sought to be minimized. In general, the 
input set D' which will minimize the LO uncertainty will be a function of the predicted 
observable. The accuracy of predictions for the triton charge radius, for example, did not 
increase (see solid and dashed red lines in fig. when B{d) instead of at was included in 
T?'. We will, nevertheless, translate the approach to the larger systems (ch. |vm| fc[D^ 
in which ^He will play the role of the deuteron as the dominant substructure in the ®He 
ground state, and as the target in elastic a-n scattering. Explicitly, LO results for the ®He 
bound-state energy and ®He properties are expected to be less sensitive to A and Arqm if 
B{a, exp) is used instead of B{t) to fit the LO three-body parameter A 3 N 1 . By reducing the 
sensitivity to Argm, I expect to obtain converged results, with respect to the basis size, at 
lower dimensionality and by that increase the flexibility of the RGM in A > 4 calculations. 
Although, the focus is on a real bound state in the six-nucleon and a virtual one in the 
five-nucleon system, even unimproved convergence will contribute to the understanding of 
the different behavior of an EET for bound- and scattering systems with regard to the 
position of the bound state in the spectrum relative to the nearest breakup threshold. To 
be specific, it is of interest to investigate what effect the relatively large gap between B{d) 
and B{t) relative to that between R(®He) and B{a) has on EET convergence. 


VII. A=4 

An investigation of systems of V -|- 1 particles with an interaction derived in the N' < N 
body sector is presented here for N = 3. The foregoing analyses of the three-nucleon 
system with EET(^) conhrmed at NLO that neither triton nor low-energy neutron-deutron 
scattering observables are universal consequences of the large two-nucleon scattering lengths. 
By increasing the number of particles by one proton, the analogous problem asks: Is the 
4-helium ground state and/or low-energy four-body two-fragment scattering observables 
on the ground states of the subsystems, i.e., deuteron, triton, and 3-hehum, universal for 
interactions parameterized by, e.g., as^t,B{t)7 
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A. Tjon correlation 


First, in analogy to the triton, correlations amongst bound-state observables are inves¬ 
tigated, specifically, the Tjon line^^ which relates the ground-state energies of the triton 
and 4-helium. By doing that, it implies the existence of a four-fermion bound state for 
all short-range interactions parameterized by two-fermion scattering lengths and the 
smallest three-fermion binding energy. While this is a non-trivial consequence itself, the 
different characeteristics of this stated compared to those of its bound two- and three-body 
fragments make it even more interesting. Specifically, it is bound by ~ 20 MeV below the 
lowest breakup threshold, triton-proton, compared to only ~ 6 MeV for the triton below 
d-n, and ~ 2 MeV in the deuteron relative to n-p. Previous LO EFT(^) calculations ^2 
suggested that this large gap does indeed follow universally with a relatively small uncer¬ 
tainty, if V = {as, at/B{d), B{t)} is used to fix the relative positions of B{d) and B{t). 
The assessment of the uncertainty was limited there to a change in input data, changing 
at to B{d), which, in the light of the results shown in fig. |^&:[^ could underestimate the 
error. In this work, we will show that not only the bound structure, but also the large gap 
to threshold can be explained, counterintuitively as it suggests large typical momenta, with 
EFT(it). Eor this purpose, I extend the LO analysisE2 to NLO in the EET(^) framework, 
and include the Coulomb interaction in the potential. The variational basis used for the 
LO and (p,n)NLO calculations includes all two-fragment groupings with bound fragments: 

G' = {triton-proton (t-p) , 3-helium-neutron (^He-n) , deuteron-deuteron (d-d)} . (59) 

Singlet-deuteron-singlet-deuteron (d-d) and di-neutron-di-proton (nn-pp) fragments are not 
included in the basis as they were found to contribute insignificantly to the bound state 
in contrast to their role in the scattering basis (see below). For those unphysical objects, 
the isospin invariance of the two-nucleon EFT(|f) up to NLO implies identical Hamiltonian 
eigenvalues for nn and d states. Both states have a total spin s = 0, an orbital angular 
momentum of / = 0, and an isospin of t = 1. As the total angular momentum and parity 
of the fragments in a certain grouping is set, it suffices to specify the maximal orbital 
angular momentum on each coordinate to define the spin and orbital angular momentum 
structure of the model space if all coupling schemes to a total angular momentum state 
= 0^ are allowed. For the LO and pNLO basis, S-waves were used on all coordinates 
except the fragment-relative coordinate between two deuterons on which, in addtition, a 
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TABLE IX: 4-helium Gaufiian width parameters 7 „ for the RGM basis in LO and pNLO calculations. 
The 7 subscripts label the Jacobi coordinate s with sVei = S 3 , and S 12 = s 1 + S 2 is the intermediate 
spin coupled from particles 1 and 2 - a neutron and a proton - within the triton state. The 81 triton, 
81 ^He, and 41 d-d states form a 1206-dimensional basis with the 6 relative width parameters. The 
width sets were derived from W 120 and optimized with the genetic algorithm (ch. E for a converged 
B{d,t) in this smaller model space. 


feG' 


widths [fm 


dim 


t-p 

& 

^He-n 


Sl2 = 1 : 71 (j) e Wg = • 




26.30585735, 

49.6384602, 26.3058574,11.7732549, 

> , 72(1) G • 

4.47594489, 

4.47594489,1.65292568, 0.81911105, 

0.81911105, 


0.33017059, 0.09380507, 0.02205199 





0.093805066 


36 


d-d 


49.6384602,11.7732549,1.65292568, 
0.330170586,0.022051993 
42.88895245,15.27652652,4.040844294, 

0.790390928, 0.138826844, 0.026116458 


S 12 = 0 : 7i(j) G Wg , 72O') G < 


71,20') G < 


45 


21 


73O) G W120, elements 3, 5, 7, 9, 11, 13 are used V/ 


D-wave was allowed, Oei(d-d) = 0, 2. This limitation in I is motivated by the analogy to 
a rectangular-well potential. Analytic bounds for the formation of a bound state on the 
depth Clo and width ^ of the well can be derivedESl which would require, e.g., for I = 1 : 
C'lot^ — 1^' Analogously, for the GauCian EFT(^) potentials with fixed LECs and 
cutoff, i.e., depth and width, particles in partial waves with higher angular momentum 
will eventually not bind together. Hence, such states will not overlap with a four-body 
bound state. Furthermore, the strength of the interaction in higher partial waves is by 
construction small in the EFT framework as shown in fig. We verified the assertion by 
comparing results for B{a) obtained in this EFT model space and one which was used for 
AV18(+UIX) predictions. The latter also includes up to three relative D-waves coupled to 
a total I = 0. The difference was found insignificant relative to LO uncertainty from the 


interaction. The basis definition is completed by the width sets given in table IX The 
nNLO basis is truncated by < 2 Vn, and the width sets resemble those defined for the 
RGM AV18(+UIX) calculationsSS] regarding their dimensionality with values optimized by 
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the genetic algorithm (app. 0 for each cutoff. 

To assess if the reduction in the number of triton components in the 4-hehum basis 


relative to its size in a three-body basis (ch. VI) is still guaranteed for an almost complete 
4-helium model space, the binding energy B{a) thus obtained was compared to one cal¬ 
culated in the much larger scattering model space. The two values differed by less than 
50 keV. This model space, used for the O'*'—channel scattering calculation, is defined below 
and is spanned by more than 7000 basis vectors. This study of reducing the dimension 
of the 3-nucleon fragments in the a particle without deviating significantly from the as¬ 
sumed converged value is crucial for applications of the method to ^ > 4 systems regarding 
practicality in terms of computing cost and accuracy in predicted observables. 

In fig. the results are compared to a LO band calculated by solving Fad- 
deev/Yakubovski integral equations quasi-exactly by |64j (referred to as FY results). The 
spread of the NLO values is not in conflict with EFT(^) which allows a < 10% uncertainty 
at NLO. Like for the triton charge radius, the upper (lower) boundary of the FY results was 
obtained by choosing different NN observables, as^t{B{d), as), to fit the LO LECs for cutoffs 
high enough so that B{a) did not change when further changing the cutoff. Figures & 


10 confirm that the LO accuracy would be overestimated with this method. The ROM 


LO calculation uses the same interactions as in sect. VI with cutoff values from 400 MeV 


(lowest red dotted line setting the lower bound for the ROM LO band in fig. 10) to 1.6 GeV 
(top solid red line setting the upper bound). In contrast to the RGM LO correlation be¬ 
tween B{t) and the triton charge radius, the width of the RGM LO Tjon band is however 
not converged. The shifts between the positions of the correlation lines corresponding to 
cutoffs of 400 MeV, 700 MeV (middle red solid line), and 1.6 GeV indicate that a variation 
of A even beyond 1.6 GeV would be necessary to assess the LO uncertainty from cutoff 
variations only. However, elaborate technical modihcations are required for the RGM cal¬ 
culation at higher cutoffs. Analogous to the conservative estimate for the LO uncertainty 
in fig. the LO Tjon correlation band is thus mapped out by both the RGM and the 
Faddeev/Yakubovsky (thick dashed lines) results, which overlap nicely in this plot. This 
combined correlation band includes the datum and the narrower NLO band, which, from 


fig. 10, has at the experimental B{t) a width of about 5 MeV centered around 28 MeV and 
results in a prediction of 


H(a,NLO) = (28±2.5) MeV , 


(60) 
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which is consistent with the expected NLO uncertainty of about 10% and with experiment. 
Comparing the NLO uncertainty to LO accuracy, which we estimate by the extremal pre¬ 
dictions at B{t,exp) for B{a), i.e. LO A = 400 MeV (lowest red dotted line) and the upper 
edge of the NLO band, to be about ±6 MeV, would be consistent with an expansion param¬ 
eter of the EFT of the same order as the naive estimate of ~ Again, the results of the 
AV18(+UIX) models he within the proposed band as it is expected of all interaction models 
of at least NLO. The observed broadening of the correlation band is a manifestation of 
the momentum dependence of the EFT expansion. The accuracy decreases with increasing 
typical momentum, eventually leading to a breakdown of the expansion. 

From the fact that there is still a one-parameter correlation, it is concluded that no four- 
nucleon contact interaction is required to renormalize the theory at NLO. One three-body 
parameter fitted to data suffices to yield proper NLO predictions for four-body observables 


within the theoretical accuracy (empty squares in fig. 10). 


B. Neutron-3-helium scattering 


As in the A = 3 system, a low-energy scattering observable should solidify the assertions 
about the power counting in A = 4 which were based in the preceding section on bound- 
state properties. In principle, all low-energy observables should be correlated with the triton 
binding energy. The recent results!^ for the singlet and triplet n—scattering lengths 
using three potential models and a N^LO chiral potential are evidence for this assertion in 
the four-nucleon scattering system. Here, the real part of the S-wave spin singlet scattering 
length ao(^He-n) for elastic 3-helium-neutron scattering is investigated. In fig. 


11 


its value 


is shown as a function of B{t) for six potentials including the Coulomb interaction. For 
two NLO potentials with A = 440 MeV (solid line) and A = 550 MeV(dashed line), with 
two-nucleon LECs fixed, we also show the effect of a smooth variation of the three-nucleon 
interaction parameter. 

To extract ao(^He-n), six two-fragment channels, ^He-n, t-p, d-d (Zrel = 0,2), (nn)-(pp), 
and d-d are included. The latter two consist of unbound fragments and usually model 
possible three- and four-body breakup reactions. They, as well as the two d-d channels, are 
however for this calculation only needed to provide configurations for distortion channels 
since only the t-p channel is open a few eV above the ^He-n threshold. For the fragment 


wave functions (see eq. (34)), a 224 dimensional basis was used for the triton and 
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FIG. 10: Correlation between the triton and ^He binding energies (Tjon line). The filled squares 
are the results using NLO EFT(]t) NN potentials from table [V| with no 3NI, and the empty squares 
represent the predictions for the 3NI fitted to B{t). The shaded NLO band results from a variation 
of the 3NI for those potentials. We compare to LO results from [^, where the upper (lower) 
dashed line was calculated with LEG fitted to Og,* Og), and to our LO RGM calculation with 

potentials from table III B{a, exp) is taken from [5S], the values for AV18 (+UIX) are reported in |49) . 
and those using a chiral expansion to NLO and N^LO in m- 


3-helium, and a 9 dimensional one for the deuteron, whose six In = 0 vectors built the nn, 
pp, and d states. For these six channels, the 20 width parameters wi 2 were used for the 


7,1 (j) in eq. (34). Almost all configurations included to build those physical channels could 
be recycled as distortion channels to allow for more freedom in the minimization of the 
variational functional. Less than ten configurations had to be excluded to avoid numerical 
linear dependences. In each distortion channel, four to six relative width parameters 'Jnij), 
taken from wi 2 with 7n(j) > 0.02 fm“^, were used. Numerical stability and convergence 
of ao(^He-n) were assessed by increasing the number of included relative widths 7^0) by 
one for each distortion channel, yielding changes in ao(^He-n) of the order of the numerical 
uncertainties, given that the initial 7n(j) were chosen appropriately. The lowest eigenvalue 
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B(t) [MeV] 


FIG. 11: The correlation between the triton binding energy and the real part of the spin singlet, 
S-wave scattering length ao(^He-n) for elastic ^He-n scattering. Blue squares represent the values 
of the NLO EFT(y-) potentials with a numerical uncertainty of 0(0.02 fm). A variation of the 3NI 
yields the continuous correlation lines for two selected NLO potentials with a fixed cutoff. 


of the Hamiltonian is equal to the ground-state energy of 4-helium in this model space and 
was allowed to change in this process by not more than 10 keV. Significantly larger changes 
in this eigenvalue which lead to a result not of the order of magnitude suggested by the LO 
Tjon band signal numerical linear dependences. If the model space is too small, or if the 
width parameters for the relative wave function were chosen inappropriately, changes of the 
order of 100 keV up to a few MeV are expected. The 3-helium-neutron scattering length 
ao(^He-n) was calculated from the corresponding diagonal element S 22 of the S-matrix at a 
center of mass energy of Ecm < 10 eV above the ^He-n thresholdl^: 


®o(-^cm) — 


1 - 522(^cm) 


(61) 


'^'\J f llf-E'cm(l + *S'22(.£'cm)) 

A fraction of the flux is diverted into the open triton-proton channel, resulting in a 
nonzero S 12 S-matrix element and hence a nonzero imaginary part of ao(^He-n). The 3NI 
variation yields an almost linearly increasing imaginary part with increasing B{t) from 
about Im{ao(^He-n)} « —7.5 fm at B{t) ^ 5.5 MeV to Im{ao(^He-n)} ss —2.0 fm at 
B{t) « 9.1 MeV. This qualitative observation is consistent with the fact that the splitting of 


trinucleon binding energies decreases with decreasing triton binding energy (see sect. VIB), 
which results in a smaller separation between the respective thresholds in four-nucleon 
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scattering. This handle on the threshold separation snggests an approach to circumvent 
numerical problems associated with the proximity of thresholds by extrapolating results for 
their physical values from calculations performed at more deeply bound 3-nucleon states. 
This flexibility is conjectured to facilitate future calculations. 

the predictions of the NLO potentials with and without 3NI for Re{ao(^He- n)} 


In fig. 


11 


decrease with increasing triton binding energy. They map out a band which includes the 
datum. A comparison to hard sphere scattering qualitatively explains this behavior. A 
higher 3-helium binding energy corresponds to a smaller nucleus analogous to the 

triton as shown in fig. As the scattering length is proportional to the radius of the 
hard sphere, ao(^He-n) is expected to decrease for increasing il(^He). We define the NLO 
correlation band to be centered around the values predicted by the potentials with zero 
3NI. The change in ao(^He-n) observed between and serves as an error estimate for 
the band centered around 7.5 fm. Again, the assumption of the center of the band at the 
experimental B{t) follows from the values in fig. 


11 


Thus, EFT(^) reports at NLO: 


Re{ao(^He-n)} = (7.5 ± 0.6) fm 


(62) 


The 3NI was varied only over the depicted range of trition binding energies, where one and 
the same fixed ROM variational space can be used. The fragment model space for the two 
interactions was optimized for the 3NI parameter fitted to B{t). Only the number of in¬ 
cluded distortion channels was changed to reach convergence within this space for each value 
of the 3NL The band mapped out by the two 3NI lines for (7.5 MeV ^ B{t) < 9 MeV) in¬ 
cludes the datum, and its slope is consistent with the one indicated by the 6 NLO potentials 
without 3NL This explicitly demonstrates that a variation of the 3NI has the same effect as 
varying the short-distance part of the NN interaction. The NN model AV18 yields a value 
within the error band. Adding the UIX three-body interaction moves this point into the 
10% NLO uncertainty radius around the datum. 

To gauge the size of the error introduced by numerical uncertainties from diagonalizing 
the Hamiltonian for the calculations with zero 3NI, ao(^He-n) was calculated over a wider 
range of energies, Ecm < 2 keV. Figure [l^ shows Re{ao(^He-n)} as a function of the matching 
energy at which it is calculated for the potential The fluctuations resulting in the 

band-like shape with a width of about 0.015 fm out of ao(^He-n) ~ 9.415 fm are due to 
uncertainties associated with the diagonalization of numerically singular matrices. The 
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FIG. 12: The real part of the spin singlet, S-wave scattering length for elastic ^He-n scattering as 
a function of the center of mass energy above the ^He-n threshold at which ao(^He-n) is calculated 


from eq. (61). Notice the scale on the vertical axis. Taking the central value of the band as the 


actual value, its decrease from 200 eV to 2 keV is due to effective range corrections, while the band is 
a result of numerical noise. was used for this plot, with the other potentials showing similar 
behavior with respect to the size of the higher order corrections and numerical fluctuations. 

error due to effective range corrections can be estimated by observing the difference, 

|Re{ao(200 eV)} — Re{ao(2 keV)}| ss 0.01 fm . (63) 

We therefore conclude that the numerical precision of our result is about 1% and is con¬ 
siderably smaller than the theoretical accuracy of ~ 10% attributed to a NLO calculation. 
For AV18, the same analysis produced stable results with respect to numerical fluctuations, 


i.e. no visible band as in fig. 1^ is found. This is attributed to the model space used for 
AVI8, which was optimized for this potential and purged of states with large but mutually 
canceling overlap with the ground state. As mentioned at the beginning of this paragraph, 
this analysis applies to values obtained in an optimized model space only. 

The conclusion drawn from the newly found correlation in fig. 0 is that every potential 
with the correct NN low-energy phase shifts and appropriately tuned three-body interaction, 
e.g., to give the correct triton binding energy, predicts not only the correct B{a) but also the 
experimental ao(^He-n) within a NLO error range. In table the results for ao(^He-n) of 
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TABLE X: RGM predictions for the triton binding energy and the spin singlet, S-wave scattering 
length ao(^He-n) for elastic ^He-n scattering of phenomenological nucleon potential models. The 
values for AV18 and AV18+UIX are taken from [55], while the AV18+IL2 numbers are new results of 
this work. The imaginary part was investigated on a qualitative level, only. 


force 

B{t) [MeV] 

Rejao} [fm] 

Im{ao} [fm] 

EFT(/) NLO 

8.48 (input) 

7.5(6) 

-2.6(?) 

AV18 

7.57 

7.81 

-4.96 

AV18+UIX 

8.43 

7.62 

-4.07 

AV18+IL2 

8.48 

7.63 

-4.28 

exp 

8.48 

7.456(20)1^ 

7.370(58)Ei] 



RGM calculations with the EFT(^) NLO potential and the AV18 NN force model with and 
without the three-body interaction models UrbanaO (UIX) and Ilhnois2 (IL2)® are given. 
The values for the recent IL2 3NI have been calculated for this work, employing methods 
described in |39l [69|. The IL2 prediction for Re{ao(^He-n)} is almost identical to the UIX 


value and not plotted separately in fig. 11 UIX and IL2 have parameters htted to, amongst 
others, the triton and 4-helium binding energies. The observed deviation of the prediction 
of both models for ao(^He-n) is therefore not a result of a dehciency of the structure of the 
potential. The correlation between B{t) and ao(^He-n) supports the conjecture that this 
small deviation can only be improved by the inclusion of higher order interactions. 


Analogously to UIX, the IL2 model easily satisfies the criteria mentioned in sect. VIA 


namely to reproduce low-energy NN observables and the triton binding energy at least with 
a 10% accuracy, and therefore its prediction for ao(^He-n) should also lie within a 10% radius 


of the datum in fig. 11 This EFT prediction is confirmed here by explicit calculation. 

Without a calculation for the scattering length in the triplet channel, the results can¬ 
not be compared to the most recent measurements^^ . The older experimental values for 
ao(^He-n) (see tab. 0 are included in the predicted universality band, and therefore we 
cannot resolve the discrepancy between the two measurements. 


60 












C. Conclusions from A=4 


The implications of these results for the scaling of four-nucleon interactions (4NI) in 
EFT(y') at LO and NLO are summarized in this section. 

Simplistic dimensional analysis suggests that a momentum-independent 4NI enters at 
N^LO, but the unusual renormalization of the 3NI may also promote the 4NI to contribute 
at lower orders to ensure renormalizability. In the RGM approach with a Gaufiian regulator, 
the cutoff was varied for the 4-helium binding energy and for the ^He-n scattering length 
from 400 MeV to 700 MeV. This led to a correlation band which is nearly identical to the 
ones mapped out when using different NN potentials (in part also at different cutoffs) or 
different 3NIs. All results were converged numerically. Moreover, we demonstrated that the 
theoretical uncertainty for B{a) decreases from LO to NLO by a factor consistent with the 
a priori expansion parameter estimate Q ~ In both cases, the physical datum lies well 
inside the NLO correlation band. 

Assuming that a 4NI enters at NLO, not accounting for it would then have two effects. 
First, the results could be unstable against cutoff variations, indicating that a 4NI is nec¬ 
essary at NLO to renormalize FFT(y'). No such effect is seen here, but a cutoff variation 
beyond the window chosen in this exploratory study needs to be performed. Secondly, the 
theoretical accuracy of the results would be reduced to LO even if results independent of 
the cutoff can be achieved, i.e., even if the cutoff can be removed to infinity without a 
4NL In other words, when one separately or collectively varies the cutoff or the NN poten¬ 
tial used or the 3NI strength, the different results should spread in a corridor set by the 
size of LO corrections, i.e., < 30%, and not by the corridor of < 10% expected in a NLO 
calculation. Above, the residual short-distance dependence of bound- and scattering-state 
observables is < 10% of the central value, taking as a conservative estimate the combina¬ 
tion of errors which occur when varying FFT(^) at unphysically short distances: Different, 
phase-equivalent NN potentials; different 3NIs; different cutoffs. The correlation bands and 
error estimates are thus quantitatively consistent with those of a NLO calculation. A resid¬ 
ual 4NI which is not necessary for renormalization but with unnaturally large coefficient 
can therefore be ruled out. 

My conclusion is that there is strong evidence that 4NIs do not enter at either LO or 
NLO in FFT(^). 
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VIII. A=5 


In the preceding section, it was shown that EFT(^) describes low-energy bound- and 
scattering- three- and four-nucleon systems systematically. In this chapter, its application 
to 5-helium is considered. The five-nucleon system is found different from its t,^He, and ^He 
predecessors. All its isotopes, i.e., ^n, ^He, ®Li, and ®Be are unstable. In this respect, 
it has no analogy in gravitational or bosonic systems. Given reasonable initial conditions, 
a system of five planets will be bound by a classical Newtonian interaction. For identical 
bosons interacting via short-range forces, bound states have been predicted for all systems 
of less than 15 particles. In ^He, both types of an interaction are present: a long-range 
Coulomb force and a short-range NN force with unnaturally large scattering lengths. The 
short-range repulsion following from the Pauli exclusion principle is implicit in the anti¬ 
symmetrized states. It is the qualitative difference in the spectra of these systems - bound 
compared to unbound - emerging from interactions with very similar characteristics that 
motivates this analysis. Besides that, it is of interest to test whether EFT(^) can be used 
for a description of low-energy five-nucleon scattering processes, a non-trivial generalization 
like the ones from two- to three-, and from three- to four-body systems. Adding a third 
nucleon made a modification of the naive power counting necessary, and there is no a priori 
argument known why this cannot happen again when a fifth nucleon is added. Yet, a five- 
nucleon force would have to be momentum dependent which suggests its promotion even 
more unlikely. The spectrum of ^He in particular contains two resonant statesE3 in the 

= I ,2 channels, at energies of ~ 0.89 MeV and ^1-5 MeV above the ^He-n thresh¬ 
old. Both energies are small compared to B(a), and the typical momenta of the neutron 
relative to the a core are then expected to be small as well. A breakdown of EFT(^) is 
therefore not expected for five nucleons. 

The universality of no shallow bound state is a different aspect. Intuitively, this non¬ 
existence of a real bound state of five nucleons is not expected to be a prediction at LO 


or even NLO in the pionless EFT. Considering the LO uncertainty in B(a) (see ch. VII) 
of about ±7 MeV, a theoretical error in i?(®He) of similar magnitude would naively be 
expected to allow for a real instead of a virtual bound state. I emphasize, this does not 
mean that 5-helium is outside the range of applicability of the pionless theory. The crite¬ 
rion of the validity of the theory is rather the prediction of a state at LO whose energy is 
renormalization-scheme independent within reasonable limits, e.g., roughly the ones found 
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in B{a). Furthermore, with NLO corrections the magnitude of this dependence must de¬ 
crease for the theory to be useful. Therefore, I search for a state of energy of about ±7 MeV 
relative to the a-n threshold. Different models of the short-distance structure of the NN 
interaction could then shift the physical two negative parity resonances below threshold to 
become bound states. If such a state is shallow, the system would be similar to the deuteron 
but contrast the three- and four-nucleon bound systems since its bound state is very shallow 
relative to the lowest breakup threshold. The (un)bound ^He would thus require a partic¬ 
ular short-distance structure and would not be a universal consequence of an unnaturally 
large pair of NN scattering lengths and B{t). 

The analysis is not intended to rival the most recent, accurate descriptions of a-n scat¬ 
tering, treated as a five-body problem. For the interaction models AV18+UIX/IL2, obtained 
using a combined variational- and Green’s function Monte Carlo technique, predictions are 
in good agreement with data for phase shifts and the resonance positions^^l. 

A comparable agreement is achieved with the combined NCSM/RGM method for chiral 
N^LO interactions^SI. The interaction in both analyses reproduced the binding energies 
of the deuteron-, triton-, and subsystems precisely, in contrast to our LO investi¬ 

gation. The essence of this section will be twofold: First, in contrast to the n-d system, 
where an accurately reproduced deuteron target resulted in an improved accuracy for the 
phase shifts (see discussion of fig. [^, a well described ^He at LO does not yield a similar 
accuracy for the a-n phase shifts. We also showed that the straightforward generalization 
of matching to a Z-parameterized amplitude instead of the less appropriate ERE does not 
improve convergence in a-n phase shifts at LO in EFT(]t). Second, the results suggest that, 
counterintuitively, an unbound five-nucleon system is a universal feature of the two- and 
three-nucleon interaction. 

The phenomena I will describe are similar to the findings for few-boson systems where 
two-boson interactions with a scattering length |a| —)• oo and repulsive three-boson interac¬ 
tions were considered^. If the triplet and singlet np scattering lengths 0^,4 are treated as 
infinite, the difference to the nuclear systems, considered here, lies in the behavior of the 
interacting particles under permutations, namely identical bosons compared to two species 
of fermions. In m one five-body bound state with a binding energy of about ten times that 
of the associated trimer was found. In addition, either a second relatively weakly bound 
state, or a state “on the verge of being bound”, was discovered. This finding is a guideline 
for our investigation. If a deeply bound state is found in the five-nucleon system, it would 
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certainly not be describable with EFT(^). The sought-for useful description of the first 
excited, shallow state with EFT(^) should not be affected by the possible existence of the 
deep ground state. Again, gravitational systems provide a heuristic justification for this 
assertion. The possible formation of a black hole from five massive objects does not render 
the Newtonian approach to the system under “normal” conditions useless. 

The analysis is comprised of two parts. First, EFT(it) interactions are used at LO to 
predict low-energy a-n phase shifts. In this course, both of the aforementioned goals are 
addressed: Is the LO uncertainty reduced if the a target is reproduced accurately? In 


this respect, does the system behave similar to d-n scattering (see ch. VIC)? Furthermore, 
the behavior of the phase shifts should already indicate the existence of bound structures. 
Second, a LO and NLO EFT(^) analysis in the ROM bound-state formalism is conducted 
to solidify the implications of the P-wave phase shift results. 

Preceding the presentation of the results, a brief, qualitative comparison of the a-n and 

the n-p system shall later allow for a better interpretation of them. One can think of elastic 

a-n, d-n, and n-p scattering as two-body problems with the distinguishing property for the 

expected behavior of the phase shifts in a given channel being the relative position of the 

thresholds and bound states. There is an important difference between the a-n and the NN 

system. Treating the a-n /j-ei = 1 amplitudes analogous to the np system would suggest 

them as higher-order quantities not subject to constraints at LO. Only the '^Si amplitude 

2 

would then be expected to be predicted within a defined uncertainty range. The shallow 
virtual P-wave state, however, discriminates the a-n from np. In [77! it was shown how the 
power counting in an a-n halo EFT for a shallow, real, or virtual bound P-wave state does 
not follow naive dimensional analysis. The halo EFT considers the a as an inert DoF, and 
hence, a-n as a two-body problem. The enhancement of a certain class of diagrams then 
parallels that in the np system in EFT(^), where it allows for the real (virtual) state in the 
(^So) channel (see ch. 0 but not in the P-waves. This similarity of the ^Pj channels in 
a-n and channels in np motivates the attempt to describe the low-lying states with 


the LO EFT potentials derived in ch. Ill for the reasons outlined below: 

First, possible bound states are expected to be formed by states with typical momenta 
between the nucleons in the core and the neutron that are small compared to momenta 
conjugate to the range of the interaction between the a and the neutron. A halo structure 
similar to ®He is expected for a weakly bound ®He, which implies a broad coordinate-space 
and correspondingly narrow momentum-space wave function. 
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Second, ^Pj bound states are the first candidates for universal properties associated 
with and B(t) or B(a), analogous to the aforementioned bosonic five-body states which 
contain relative S-waves only. Universal shallow ^He states cannot be pure S-waves, because 
such states are not accessible for three identical fermions. As I will argue below, the a core 
is described in terms of an S-wave basis at LO and pNLO. Therefore, one must allow at least 
one relative P-wave for a totally antisymmetric wave function. Furthermore, the P-wave 
channels are considered here as a subsystem of the ®He state with its two halo neutrons in 


relative P-wave states but a total angular momentum of zero (see ch. IX). 

The construction of the RGM model space is presented in some detail here, because of its 
peculiar simplicity compared to the larger spaces employed in potential-model calculations. 
Already in the deuteron, the deviation from most approximations to its wave function is 
apparent since it is a pure S-wave state up to NLO in EFT(y'). The construction of the 
^He RGM model space starts with selecting an appropriate set of groupings of the two 
protons and three neutrons. A a-n division only is used here. As we consider the a-n 
scattering system only for Ecm < 10 MeV, i.e., below the lowest breakup threshold, this is 
a reasonable approximation. For the hypothetical bound state, more care has to be taken. 
Assuming that such a state is similar in its shallowness, momentum distribution, charge- 
and matter radius to the ®He halo structure, the approximation would certainly be justified, 
too. It is important that this assumption contradicts the possibility of a similarly large LO 
uncertainty in the energy of the four- and five-nucleon system. According to this estimate, a 
bound state of about 7 MeV below threshold would be admissible but certainly not exhibit 
a pronounced halo structure. Nevertheless, we will justify a posteriori that other groupings 
like t-d are negligible. The reasoning is the following: if such a fragmentation is significant, 
so will states of the a-n structures that resemble relatively small separations between the 
core and the neutron - both objects are small. This is not found in our calculations, and 
the predictions are derived with the a-n group structure only. 

Of the two fragments, only the a has a nucleon substructure. Its four constituents 
can be divided into two deuterons, or a triton and a proton, or a ^He and a neutron, or 
unbound nn-pp and d-d spin singlets. This leaves the two- and three-body matryoshkas to 
be opened: the triton reveals a triplet and singlet np system with a neutron, and the 3- 
helium just replaces the single neutron with a proton. Thereby we have defined all relevant 
groupings and the accompanying angular momentum coupling schemes and definitions of 
Jacobi coordinates. The specific values for the orbital angular momenta and the GauBian 
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width parameters are now to be chosen to appropriately represent the small dolls in the 
larger ones, i.e. for example, the triton within the a. 

The lightest bound structure is the deuteron, which is well described as an S-wave state 
up to NLO (ch. 0 . The bound A = ?> systems are therefore in a first approximation modeled 
by two S-waves. The angular momenta on the two Jacobi coordinates of the three-body 
system have to be the same, in any case, to couple to a total I = 0. Only states with total 
I = 0 states are of interest, because only they have nonzero overlap with the four-nucleon 
ground state in this approach. From a technical perspective, this is because the interaction 
is even at pNLO a spherical tensor or rank zero, i.e, a central operator that does not mix 
states with different total orbital angular momenta 1. Contributions to the bound three- 
nucleon binding energies at LO EFT(y') from states with two P-waves were found to be less 
than ~ 100 keV. Even relative to NLO accuracy, this is a small amount, and the omission 
of such configurations is justified. Now, with d, t, and ^He as S-wave states, an S-wave for 
the third coordinate of the enclosing “^He doll is implied. Using a short-hand notation, the 
ansatz for a ^He wave function used in this calculation is: 

I^He) = |d-d)ii8|n)i4 -b |t-p)85|n)i4 + |^He-n)87|n)i4 . (64) 

The subscripts denote the number of basis vectors. For instance, there are 
85(triton-p) x 14(relative a-p widths) t-p configurations in a five-nucleon basis basis of 
dimension 4060. Fourteen different width parameters, 

f 1.6902,1.3423,0.8214,0.4447,0.2574,0.1690,0.1185 1 

W„_n = < > 

[ 0.0843,0.0500,0.0257,0.0186,0.0097,0.0056,0.0028 J 

expand the relative motion between the neutron and the 290-dimensional a. With the 
20-dimensional deuteron wl20 as base, and the same 20 width parameters to expand the 
d-d, d-n, d-p, d-p, and ^He-n functions, we selected two sets of 290 basis vectors from the 
8000 for ^He. One set was optimized to yield maximal B(a) for A = 200 MeV, the other 
for A = 1600 MeV. The two model spaces of equal size define the values Aj^gm,{i, 2 ) used 
in fig. [T^ For all A, the calculations are conducted in both model spaces as an additional 
assessment of the renormalization-scheme dependence of the phase shifts parameterized by 
these two Arqm values. 
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FIG. 13: Phase shifts for elastic a-neutron scattering. Results of LO EFT(-^) potentials with cutoffs 
from 300 MeV - upper (lower) edge of the shaded area in the left (right) panel - to 1.6 GeV - lower 
(upper) edge of the shaded area in the left (right) panel - are compared to R-matrixI^ fits. The 
shaded areas mark the assessed LO uncertainty by a variation of four renormalization parameters: 
the regulator parameters A and Argm (corresponding to two model spaces for the a core, black and 
red lines), the three-body interaction strength A 3 N 1 , and two fit input sets (solid lines use B{t), the 
dashed line uses B{a) with all other renormalization parameters fixed). 


With the variational basis defined, I will discuss the LO treatment of elastic a-n scat' 
tering in the and ^Pj channels. The leading-order potentials specified in table 


III 


are 


employed, thereby covering a cutoff range from 300 MeV to 1.6 GeV. In order to search 
for bound states, I am looking for a rapid decrease of the phase shift in the respective 
channel through J (see, e.g., 178]) and a corresponding eigenvalue of the Hamiltonian below 
the lowest, i.e., a-neutron, threshold. This threshold is identical to B{a), and varies with 
Argm, he., the model space used for the a fragment. 


The results for the phase shifts are shown in hg. 13 Within LO accuracy, the phases 
do not indicate a bound state in either channel. No crossing through ^ is observed. The 
shaded areas are mapped out by all our LO predictions and assess the uncertainty by their 
maximal deviation from the R-matrix data (blue). The spread was induced by the different 
short-distance structure of the interaction. The various models are parameterized here 
analogous to how we probed for renormalization-scheme dependence in A < 4 systems: 
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1. Variation of the regulator parameter A from 300 MeV (upper (^S'i)/lower (^Pj) edge 

of the shaded area) to 1.6 GeV (lower (^S'i)/upper C^Pj) edge of the shaded area); 

2 

2 . Change of the variational basis, i.e., Arqm, using two a core basis sets (solid 
red/black) with which B{a) is predicted consistently within LO and pNLO uncer¬ 
tainty; This corresponds to the above mentioned reduction of the basis size from 8000 
to 290 vectors, without significant effect on the binding energy. 

3. Two sets of input data to determine the LO LECs: 

= {as, B{d), B{t)} (solid lines), = {ag, B{d), B{a)} (dashed lines); 

4. Various values of the three-nucleon interaction parameter A 3 N 1 for which the resultant 
B{a) and B{t) differ by less than 30% from their experimental values. 


The third option is to be compared to the d-n system. Here, it is the a which is either 
predicted less strongly bound (V[) or taken as input (P^)- contrast to fitting to B{d) 
in the three-nucleon case, no improvement in the predicted phase shifts is found for an 


A 3 N 1 fitted to B{a,exp). On the contrary, the example displayed in the left panel of fig. 13 


shows for A = 300 MeV and within one model space that using T >2 yields a larger deviation 
(red dashed) from data relative to (red solid). This confirms the difference to n-d 
scattering where an accurately reproduced deuteron fragment improved the accuracy of 
phase shift predictions significantly. Regarding the fourth handle defined in the list above, 
we found the phases to converge to the upper edge of the shaded area when A 3 NI was varied 
such that H(q;,Argm) is increased even beyond H(Q;,exp). At a certain A3NI, the system 
develops a resonance and the graph of the phase shifts discontinuously changes its shape, 
dropping rapidly through | at about 1 MeV. This signals the appearance of a bound state 
for a sufficiently weak 3NI. I assume this state to be the analog of the deep ground state 
found for bosons, and therefore refrain from a further investigation here, because of the 
inappropriate model space and the expected breakdown of the pionless theory for such a 
deep state. 

Using B{d) instead of at for the determination of the LECs did not change a-n phase 
shifts significantly, although the respective B{a) differ by up to 5 MeV. The small difference 
in phase shift predictions for the two Argm values (solid black and red) could be compen¬ 
sated by a variation of A 3 NI; thus verifying once more the interpretation of an appropriate 
truncation of the variational basis with a regularization parameter. Eurther dependence 
on Argm was probed by changing the number of included distortion channels from 800 
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FIG. 14: Energy levels of the ^He system for = 5 "*", J calculated with the RGM for EFT(y-) 
interactions. For three values of the regulator parameter A (eq. the LO (dashed) and pNLO 

(solid) predictions are displayed. The a-n thresholds are set by B{a) and are predictions at LO 
(gray) and input at pNLO (red). 


up to 1280, with no significant effect. This justifies a posteriori the omission of the other 
five-body groupings - e.g., the closed t-d channel - in the basis which are not expected to 
contribute as distortion channels when channels which are already open yield only small 
improvements. 

The predicted phase shifts in the ^Pi and channels are identical because the LO 

2 2 

interaction is independent of J and J . The power counting suggested by halo EFT for 
shallow P-wave states includes L ■ S operators at loEI. Consequently, the two P-wave 
channels are coupled; a feature that would express itself not before N^LO in the NN EFT(^) 
considered in this work. Therefore, I do not attempt to refine predictions for scattering 
phases with a pNLO EFT(^) calculation, as it does not yet include J dependend operators, 
either. Compared to the np P-wave phase shifts (fig. [^, the a-n analogs are also small but 
are predictions and not input. 

the results of RGM bound-state calculations in the and ^Pj channels at 


In fig. 


14 


LO and pNLO in EFT(^) for three cutoff values, A = 300 MeV, 400 MeV, and 500 MeV, 
are presented. The discussion in app. [D] explains this limitation to relatively small values. 

When the model space was defined, it was already mentioned that a ‘^Si state of five 

2 

nucleons is Pauli forbidden at LO because, at this order, all nucleons are in relative S-wave 
states, i.e., the coordinate-space wave function is symmetric like the isospin function of 
three nucleons of the same species. As there is no totally antisymmetric spin state for 
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three spin-i functions, the ^Si five-nucleon state vanishes. Although this is precisely what 


we confirm for A = 300 MeV and 500 MeV in fig. 14, a detour on how such states could 

emerge as artifacts of the finite RGM basis is given. An asymmetry is introduced in the 

coordinate wave function by using different width parameter sets to expand the radial wave 

functions of the various Jacobi coordinates. This will allow a ‘^Si state, nevertheless, in 

2 


the sense that the explicit antisymmetrization of the RGM basis (eq. (42)) results in mixed 
symmetry states in coordinate- and spin space. Only if a sufficiently large basis is chosen 
for all Jacobi coordinates, a then symmetric coordinate space wave function will ensure 
that antisymmetrization will eliminate this state. In practice, this limit can be identified 
by numerical instabilities due to numerically zero eigenvalues of the norm matrix because 


of the almost zero states. This leads via eq. (44) to large eigenvalues of the Hamiltonian 
which in turn are sensitive to modifications of the basis at the level of adding/removing 
single vectors. In effect, the asymmetry introduced by the small basis models, e.g., states 
where the a fragment includes components in its basis with two I = 1 angular momenta, 
still forming a state of total I = 0. This suggests an interpretation of a converged B{a) 
in a relatively small RGM basis. The explicit addition of structures with two relative P- 
waves could have been found small, partly, because they are simulated by an appropriately 
distorted pure S-wave basis. The four-body calculation was carried out in a relatively 
large model space, in which a totally symmetric coordinate space wave function can be 
assumed. Reducing the number of basis states by removing those which do not decrease 
B(a) signihcantly, introduces the mixed symmetry states in the ^He basis. The important 
conclusion from this is not to exclude basis states, which an EFT power counting would 
demote to higher orders, a priori, when operating in small RGM model spaces. Therefore, 
I will include conhgurations in ®He with I = 1 on the coordinates of the two (halo) neutrons 
coupled to a total I = 0. An incomplete basis with all five relative coordinates in S-waves 
introduces such configurations implicitly but will also yield a bound system. Therefore, 
bases with both P- and S-wave vectors bear more potential for the formation of numerical 
linear dependencies if the S-wave configurations with more distorted S-wave configurations. 
The extreme scenario would model the six-nucleon state with a distorted pure S-wave basis 
only. This allows a severe reduction in computing time while exploiting the insignificant 
contribution at (N)LO for nucleons in relative P-waves utilizing a property of our EFT(^) 
potentials. 


70 





Returning to the discussion of fig. 14, the -calculation yields for both cutoffs an 


eigenvalue (dashed (LO) and solid (NLO) black) slightly above the a-n threshold - deter¬ 
mined by B{a) at the respective order (gray (LO) and red (NLO)). The corresponding state 
is identified as an a. core with an unbound neutron. This interpretation is verified by the 
disappearance of the state from the spectrum, once the broad basis states which expand the 
continuum relative wave function between the a and the neutron are omitted. The corre¬ 
sponding eigenvalue is not to be confused with the position of a potential hve-nucleon a-n 
resonant state. Information about the latter mandates boundary conditions on the relative 
wave function between the “^He and the neutron appropriate for a scattering state. Such a 
condition (eq. [M] ) is not enforced here, and although the state might have an overlap with a 
resonant state, no information about the latter can be deduced. For that purpose the RGM 
scattering formalism will be used in a future analysis. The bound state pNLO calculation 
in this channel is then just another, yet powerful demonstration of the Argm independence 
of the pNLO potentials as B{a, exp) is accurately reproduced here within a different model 
space than the one used to fit the LECs. 

While this is the result expected from the Pauli principle, the qualitatively identical 
results in the ^Pj channel are non-trivial with a profound implication. For all three cutoff 
values the spectrum does not contain a bound state even in this channel. This finding sug¬ 
gests that the unbound nature of a five-nucleon state does follow universally from unnatural 
and B{t). Explicitly, I conclude from the results, that a system of hve fermions does 
not sustain a shallow bound state when all particles interact via short-range forces charac¬ 
terized by and B{t) and, additionally, two particles repel each other with a Coulomb 
force. This does not stand in contrast to the bosonic case. A deeply bound state, as it 
was found there, might also be an element of the ®He spectrum but have been cut off by 
the employed RGM regulator. In more detail, a state with a large binding energy relative 
to the a is also expected to be accordingly more localized in coordinate space. Therefore, 
its expansion in terms of Gaufiian functions requires narrow width parameters. As those 
parameters are absent in the employed basis, the small structure of the hypothetical deep 
ground state cannot be resolved. 

Finally, a limitation of the analysis presented here shall be mentioned. The universality, 
i.e., the independence with respect to details of the microscopic structure of the interaction, 
of the unbound five-nucleon system was inferred from three models of this short-distance 
structure only, corresponding to the three cutoff values. Yet, in combination with the 
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FIG. 15: Energy levels of the ^He system in a ^Pj state calculated with the RGM for EFT(^) 
interactions without (left) and with (right) Coulomb interaction. For a regulator parameter of 
A = 400 MeV (eq. ([20|), the LO (dashed) and pNLO (solid) predictions are displayed. The a-n 
thresholds are set by B{a) and are predictions at LO (gray) and input at pNLO (red). The left 
column uses the same potential as employed for fig. while the right one refits the 3NI parameter 
A 3 NI to yield i3(a,exp) even without Coulomb repulsion. 


phase shift results, even this limited probe provides strong support for our conclusion. The 
dependence of the phase shifts on the renormalization-scheme was analyzed much more 
thouroughly and results were found consistent with the nonexistence of a bound state. For 
instance, the phases predicted with A = 1.6 GeV do not differ to an extent that would 
suggest a profound difference of the spectrum relative to that of, e.g., A = 800 MeV. 

This would contrast the appearance of a bound state for cutoffs or other renormalization 
parameters for which no explicit RGM bound-state calculation was conducted. Instead of 
extending the analysis in this direction, we therefore deem the assessment of the role of 
the Goulomb repulsion - whether it is decisive for rendering the system unbound - more 
interesting. The contribution from the Goulomb interaction to the energy eigenvalue was 
throughout of the order of -|-1 MeV, i.e., the Coulombless system is more deeply bound. 
Intuitively, one could imagine that it is this small amount which makes the difference 
between the bound- and unbound system, especially because such an interaction was not 
considered in the calculations for bosons. The predictions calculated with the proton charge 


set to zero (left spectrum) are compared in fig. 15 to the previously discussed results (right 
spectrum) with charged protons. Protons and neutron are still distinguishable but the 
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interaction is now a scalar in isospin space, i.e., nn-, np-, and pp interactions are identical. 
We used the same model space and one fixed regulator with A = 400 MeV. At LO, this 


results in a four-body bound state (gray line in fig. 15) ~ 0.91 MeV below the B{a,LO) 


result of fig. 14 This energy sets the lowest breakup threshold for the five-particle system, 
whose lowest LO eigenvalue was calculated as i?(5,LO) = 23.27 MeV (dashed line). Using 
the same potential as before with nonzero Coulomb interaction, the pNLO calculation no 
longer yields B{a,exjp) but a deeper bound four-body state at about 28.88 MeV (red solid, 
left column). In effect, not the binding energy of the a was used then to fit the LEG but this 
shifted value. This is admissible because there is no intentional experimental measurement 
for an uncharged four-nucleon system. Any value of about the same distance from threshold 


is appropriate. Nevertheless, the right column in fig. 15 shows the results for a 3NI whose 
parameter was adjusted at NLO to yield a four-body binding energy of 71(0;, exp). In both 
approaches, the predicted lowest energy eigenvalue of the five-body system (solid black 


in fig. 15) was about 0.56 MeV above this threshold at pNLO. This result suggests the 
generalization of our previous conjecture: A system of two fermions of type one plus three 
fermions of type two which interacts via short-range forces characterized by two unnatural 
scattering lengths and a three body binding energy, does sustain a four-body bound state 
but no shallow five-body bound state in its vicinity. The Coulomb repulsion between two 
of the five particles is not responsible for this unexpected phenomenon. This finding is to 
be confirmed for a probe of the sensitivity of the difference |71(4) — 71(5) | towards a change 
in the short-distance behavior of the interaction. Here, the Coulombless case was only 
considered for one cutoff value and two renormalization conditions. However, compared 
to the “charged” case, the threshold separation did not change signihcantly and naively it 
would come as a surprise if this pattern would change with A or Argm- 


IX. A=6 

Within the limits regarding cutoff variation and model-space dependency, it was shown 
in the previous chapter that the physical unbound ^He system is a prediction of the pionless 
EFT. In this section the aim is to predict not only qualitative features like this but also to 
demonstrate the convergence of the EFT for the next heavier nucleus, 6-helium. The study 
is exploratory, and for definite conclusions about the usefulness of the pionless theory for 
the description of the 6-helium ground state we rely on computational resources unavailable 
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at present. Nevertheless, the results demonstrate the insight the EFT approach is able to 
provide. Furthermore, the variational RGM method is established as another numerical 
tool for the investigation of extended halo structures (see e.g. [80] for an investigation using 
the hyperspherical-harmonics and the coupled-cluster methods). 

In contrast to 5-helium, ®He sustains a real shallow bound state and not just shallow 
resonant states. Besides lying on opposite sides of the breakup threshold, the two systems 
are similar with respect to their pole position relative to the lowest breakup threshold. 
Experimentally, the J'^ = O'*' ground state of ®He is only bound by about 0.975 MeV below 
the three-fragment a-nn breakup threshold, set here by B{a). This value is of the same 
order as the pole positions of the lowest ^Pj ^He resonant states above threshold. It is small 
relative to the binding energy of the a fragment in the lowest breakup channel. 


We have shown in ch. VII that the a system can be described with EFT(y'), and because 
the two additional neutrons increase the energy in the system only marginally, typical 
momenta of nucleons in the six-body system are also expected to be below the breakdown 
scale of the pionless theory. Therefore, it is a worthwhile endeavor to apply EFT(^) to 
the six-nucleon system with the expectation that the physical state is within its range of 
applicability and hence predicted with increasing accuracy order by order. This would be 
a remarkable consequence because of the peculiar halo structure of the ®He ground state. 

The attribute “halo” refers to the existence of sufficiently distinct momentum scales in 
the system ED: 

• In an A-nucleon system, one scale (M^i ~ is associated with the excitation 

energy Ec of a core structure. This core is composed of a subgroup oi A' < A particles. 


• A second scale (Miow ~ ^'m^Eh) which is small relative to Mhi and is set by the 
separation energy Eh of nucleons (one or more) from the bound A-body system. 

As a consequence of the low separation energy, the probability to detect neutrons at dis¬ 
tances large compared to the size of the core is unusually high. The large extension of the 
halo nucleus has in turn important consequences for low-energy scattering reactions. The 
corresponding cross sections are then found very different from those where the target’s size 
is given by the interaction range of its constituents. Another peculiarity of the ®He halo 
in particular is its borromean structure, i.e., its lowest breakup channel consists of three 
fragments. That corresponds to the three interlocked rings which separate if only one of 
them is broken. Analogously, if one neutron is removed from ®He, the remaining a and 
neutron will not form a bound system. 
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According to the given definition of a halo, we already discussed two halos in this work: 
first, the deuteron as the lightest nuclear halo. Consistent with the pionless theory, the 
proton core of this two-body state can be treated as a point particle for probes which transfer 
only momenta less than Mhi ~ m-,^. Compared to this Mhi, the separation energy of the halo 
neutron, namely the deuteron binding energy B{d) ~ is certainly small. The triton also 
meets the criteria and is an example of a two-neutron halo state. Again, the proton-core 
scale is large relative to the neutron separation energy of about ~ B[t)—B{d) « 6.3 MeV. 
The a particle does not exhibit a halo structure. Its neutron separation energy of about 
21 MeV exceeds the excitation energy of the remaining 3-helium. With all five-nucleon 
systems unbound, ®He is the next largest, stable halo nucleus. Like the triton a two- 
neutron halo, its core structure is an a particle with an excitation energy of approximately 
Ec ~ 10 MeV set by the t-p threshold. Compared to this Ec, the neutron separation energy 
of 6-helium E^ ~ 1 MeV is reasonably small. 

In contrast to the deuteron and the triton, the halo character of 6-helium would be a 
prediction of EFT(^). In the two-nucleon system, a finely tuned short-range interaction 
produces unnaturally large scattering lengths and by that the halo. Without this fine 
tuning, e.g., between strength and range of the interaction, a generic potential is more 
likely to produce a scattering length of the same order as the interaction range, which 
would correspond to a deeper lying ground state. The shallow triton constrains the nuclear 
interaction further. Not all interactions that produce a shallow deuteron predict a shallow 
triton. A particular three-nucleon parameter, or a specific cutoff value, as two specific 
examples for renormalization parameters, have to be chosen for the reproduction of the 
physical, shallow triton state. The two- and three-nucleon halos are therefore input for 
EFT(|t), while the six-nucleon halo would, without modification of the power counting, be 
a prediction. The halo character might, however, express itself only at higher orders, which 
would make it a consequence of details of the interaction, e.g., specific values of higher-order 
effective-range parameters {ri and vi in eq. and not a universal property of the shallow 
deuteron and triton B{t)). This is because of the LO uncertainty assessed for B{a) 


in ch. VII of about 7 MeV. An uncertainty in B{6) of similar magnitude would allow for 
neutron separation energies of comparable size to the approximate 10 MeV Miow scale set 
by the a core. Therefore, the clear separation of scales would develop at higher orders, 
where the usefulness of EFT(^) is determined by the amount by which the uncertainties of 
predictions for \B{6) — B(a)\ decrease from LO to NLO. 
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FIG. 16: The ^He-neutron-neutron group structure of the RGM basis for 6-helium. The theoretically 
predictedl^ opening angle 9 = (83°)rt;Q° motivates the choice of width parameters 7 i,(i, 2 )- 


We employ RGM bound-state calculations in the J'^ = O'*" 6-helium channel at LO and 
pNLO within EFT(^) in order to assess whether a bound six-nucleon system is a universal 
consequence of shallow two- and three-nucleon states, and if the halo character of this six- 
body object is a perturbative effect. First, the RGM model space is defined. I follow the 
same procedure as in the previous chapter: 

• selection of relevant groupings/fragmentations G' which determine the 
Jacobi coordinates; 

• selection of maximal angular momenta for the Jacobi coordinates and accompanying 
orbital- and spin-coupling schemes; 

• selection of sets of Gaufiian width parameters appropriate for the expected spatial 
extension of the six-nucleon state in the range of applicability of EFT(^). 

Gonsistent with the physical two-neutron halo structure around an a core, previous RGM 
analysesS3 with model interactions identified the an-n fragmentation as the dominant 
structure. This ^He grouping alone was found to bind ®He. Other group structures, e.g'., 
triton-triton, typically contributed less than 1 MeV to the ground-state binding energy. An 
expected LO uncertainty of about 10 MeV justifies not to include these configurations in 
the variational basis. For this work, we therefore approximate the ®He ground state by the 


single Y-shaped ®He-n group (left sketch in fig. 16). The other Q!-2n configuration (T-shape, 


right sketch in fig. 16) differs from the Y-configuration only in the definition of the Jacobi 
coordinates. The truncation of the partial-wave expansion of the two neutron coordinates 
rei ^2 discriminates the two structures. They are equivalent if an infinite (complete) set of 
width parameters and angular momenta is used to expand the functions dependent on ni^ 2 - 
The contribution of T-groups is estimated to be of the same order of that which is due to an 
increase of the number of widths in the Y-groups. I choose the sets from which the 7i,(i,2) 
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widths (defined in fig. 16) were taken such that adding new widths to those sets did affect 
the eigenvalues of the LO Hamiltonian by less than 1 MeV. These are the widths used for 
approximating the neutron halo around the core. An explicit inclusion of T-groups is hence 
unnecessary considering the expected accuracy of a (N)LO calculation. 

The angular momentum structure is defined next. For the a core, the same arguments as 
given in the ®He chapter justify a pure S-wave ansatz of t-p, ^He-n, and d-d groups. For the 
two neutron coordinates, the respective angular momenta have to be equal because the total 
angular momentum has to be zero. Again, this is because the central (pN)LO interaction 
cannot induce transitions between states differing in their total angular momentum. Hence, 


two sets (/i ^2 = 0,1 in fig. 16) of angular momenta are considered: one with two S-waves 
and the other with two P-waves coupled to a total L = 0. 

An appropriate choice of width parameters for 7i^(i^2)) halo neutrons, and the 

a core completes the definition of the variational basis. Initial calculations were carried 


out in a relatively large model space (MS(large) in top graph of fig. 17). This basis could 
be reduced to span a smaller model space (MS(small) of dimension < 720; results in all 
graphs of fig. First, I define the large space before its reduced derivative is intro¬ 

duced. Heuristically, such a basis should be able to expand an a structure whose binding 
energy, and thereby approximate spatial extension, is known at each order of the calcula¬ 


tion. From the large four-body model space used in ch. |VH[ a subset of 85 t-p, '^He-n, and 
117 d-d components was selected such that B(a) calculated in this 203 dimensional space 
was only about 1 MeV lower than in the large space. The parameter sets of the radial 
part of that component of the wave function which depends on ni ^2 were chosen under 
the following considerations; Basis vectors which correspond to small separations between 
the core and a neutron are expected to have little overlap with a state within EFT(^)’s 
range of applicability. Such highly localized states are associated naively with large typical 
momenta. We therefore allow only width parameters up to a certain 7max as larger widths 
correspond to narrower Gaufiians. I choose the cutoff width for each regulator cutoff as: 
7max ~ 1 fm“^ -|- Starting from this 7max, three values were defined as means of normal 
distributions. The three distributions were used to generate 3 x 203 width parameters, 
where each set can be imagined as a compact a core with an attached neutron in a shell 
whose distance from the core is set by the respective mean value. For the widths of the posi¬ 
tion of the second neutron, the same procedure was implemented but with five distributions 
generating 5 x 108 different parameters. The smallest of those five central values was chosen 
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as 7min = 0.0001 fm“^. Vectors including widths from this distribution correspond to the 
broadest GauBians in the basis which are necessary to allow for the experimentally®®^ 
and theoretically® established correlation between the neutron coordinates, which places 


them at almost opposite sides of the core structure (hg. 16 right T-configuration). 

I summarize the model space definition by giving the total number of basis vectors: A 
203-dimensional a core for each of the two sets of angular momentum coupling schemes 
for the remaining two neutron coordinates; three sets of width parameters expand the a-n 
coordinate, each covering a different range between the two objects; five width sets expand 
the an-n Jacobi coordinate, allowing for both, an extended halo structure as well as a more 
deeply bound, localized state. Hence, the basis is of dimension 2 x 203 x 3 x 5 = 6060 m- 

A smaller basis was sought to facilitate the analysis of numerical instabilities associated 
with the RGM method, e.g., sensitivity of predictions to specific basis states correspond¬ 
ing to either very broad or localized wave functions, the influence of the regularization of 
irregular Coulomb/Bessel functions, or the dependence on the energy used to extract the 
scattering lengths for the LEG fit. For instance, we could sort out pNLO LEG sets which 
were numerically harmless up to the four-body problem but are unstable in the six-body 
system due to the increase of the expectation values of single operators. It was also possible 
to assess much quicker the impact of a shift of the interval from which the width parame¬ 
ters for the halo neutrons were selected. This enabled us to optimize this range such that 
the second lowest eigenstate was reproduced more accurately in the small basis. Compare 


values marked by the full/empty circles in the upper graph in fig. 17 for a dimension 288 
and 4060. The value obtained in the much smaller basis is converged to a lower value due 
to a wiser selection of the interval for 7 j^(i^ 2 )- 

The smaller model space was derived from the large, defined above, by using the same Y- 
group structure with its two angular momentum coupling schemes but reducing the number 
of vectors responsible for the expansion of the a core and the first attached neutron. This 
basis consisted only of 144 x n vectors, where n is the number of width parameters used to 
expand the dependency on the coordinate of the second neutron. The results are presented 
for n = 2,3,5 in hg. with each increment adding more localized basis states, only. 
For A = 300 MeV and B{a) G ^lO’ found that this reduction from 6060 to 720 

basis states affected the lowest eigenvalues insignihcantly or even predicted lower, better 
converged values. 

The EFT interactions we utilized in this section are specihed next. To assess the sensi- 
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FIG. 17: EFT(^) results for the (second) lowest energy eigenvalues ((circles) squares) in the 6-helium 
channel as a function of the basis size. (pN)LO predictions ((full) empty) are shown for three cutoff 
values with either B{t) (A = 200 MeV) or B{a) (A = 300,400 MeV) as experimental input to fit 
LO LEG. For A = 300 MeV, a small (left upper graph) and large (right upper graph) model space 
was analyzed including values for a particularly (broad) localized basis ((diamonds) triangles). The 
red and gray horizontal lines mark the a-nn breakup thresholds. 


tivity of predictions on particular models of the supposedly unobservable short-distance 
structure of the theory, three regulator values were chosen, A = 200,300,400 MeV. 
The reason for using relatively small values for the cutoff, were the relatively ac¬ 
curate predictions for B{a) if B{t) was LO input, e.g., for A = 200 MeV and 
B{t,LiO) = B{t,exp) B{a,hO) ~ 30 MeV. For increasing A, a smaller value for 


B{a,LO) was predicted, deviating by as much as 7 MeV from experiment (fig. 10). As a 
consequence, the necessary pNLO corrections to B{a) increase with A, while corrections 
to at, one of the other constraints used to determine the pNLO LECs, remain relatively 
small. Specifically, a 7 MeV pNLO contribution constitutes about 25 % of B{a), while 


typical corrections to at (condition (ii) in ch. Ill) do not exceed 5 % of the physical at- 
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The magnitude of the resulting LECs was found to increase rapidly with A in consequence. 
Correspondingly, the contributions from individual operators increase, too, only to cancel 
in yielding a small number. Here, I avoided those numerical difficulties by using small cutoff 
values but stress that one is not bound to them, neither by the RGM, and certainly not by 
the EFT framework. 

NLO corrections were included perturbatively with LECs determined as described in 


ch. Ill The specification of the LO fit input T>' concludes the discussion of the interac¬ 
tions. We used ~ { 0 ‘s,B{d),B{a)} to determine the LEC for A = 300,400 MeV, and 
"Dlo ~ Ws, B{d), B{t)} for A = 200 MeV. We used B{a) as input for the cutoff values 
which would have predicted a value for it deviating more than 4 MeV from experiment 
if B{t) is used instead. This was done to reduce the pNLO correction to A > 2 binding 
energies in order to yield LECs which kept individual operator matrix elements of about the 
same order as the binding energies. Thereby, the numerical calculations became feasible. 
We could use B(t) for A = 200 MeV because its prediction for B(a) was off by less than 
2 MeV. Chosen for practical reasons, this input is just another way of renormalizing the 
theory and does not violate established EFT regulations. 

Figure [T7| displays the predictions for the lowest eigenvalues of the LO Hamiltonians in 
the two-proton four-neutron O'*" channel for various basis dimensions. For A = 300 MeV 
(upper graph), the values were calculated in both the small and large model spaces defined 
above. The larger basis does not place significantly lower bounds on the binding energy. 
We even find the smaller basis more appropriate to expand the states closer to threshold 
because of the refined choice of halo-width parameters (see above). As the cutoff is only 
varied by ±100 MeV and does not exceed the constant RGM regulator for the asymptotic 
functions (eq. 39, (3 = 1.1 fm~^) significantly, the other calculations were carried out in the 


smaller model space. Increasing the dimension means adding basis vectors with larger and 


larger width parameters on the n 2 coordinate (hg. 16). By that, more and more localized 
structures are included. The stability of our results confirms that even the smallest bases 
are apt to expand the states to small enough distances. The newly added vectors are thence 
identified with such localized states, which represent short-distance structure that should 
according to EFT(^) tenet not contribute significantly to observables. As it cannot be 
assumed that a halo structure emerges in EFT(it) at all, it is imperative to increase the 
basis further to allow for more localized states, even if a certain eigenvalue is found stable. 
However, these eigenstates will eventually assume energies too large for a description with 
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EFT(^). The states close to threshold - above or below - are naively thought to be in the 
range of applicability. 

To further the understanding of the interplay between EFT and RGM regulators at LO 
and NLO, I selectively removed basis vectors that correspond to relatively short distances 
between nucleons in the a core. As expected, the LO predictions (upper graph, empty 
diamonds for dimension 3150) are almost unaffected due to the relatively long-range nature 
of the interaction at LO. At NLO, the short-distance structure becomes more relevant. From 
a technical perspective, this increasing importance resembles matrix elements of respective 
operators which receive most of their contribution from small distances. Consistently, the 
NLO correction for the shallower of the two states is found considerably different, albeit 
the dimension is practically the same compared to the prediction of a basis without this 
deficiency (compare filled diamond and circle at dimensions 3150 and 3248, respectively). 
Without the short-distance components the basis cannot expand states irrelevant at LO 
but of apparent significance at NLO. This concludes the discussion of the convergence in 
the variational basis. 

For A = 200, 300 MeV, two bound states (empty circles and squares) below the a- 
nn breakup threshold (red/gray lines) are found at LO. Referring to the case of identical 
bosons interacting via comparable, short-range forcesES, multiple bound states with binding 
energies differing by orders of magnitude do not come as a surprise. However, such states 
are expected to be outside of the range of applicability of EFT(if), i.e., exhibit a strong 
A and T)' dependence, while those states which are in the range are excited states in the 
spectrum. 

Naively, the deeply bound states should not be in the spectrum at all if the regularization 
truly eliminates modes of high momentum, as suggested. It is the variational basis through 
which such high-momentum modes are reintroduced, because I am at present not able to 
derive bounds for the Gaufiian width parameters which define the basis from A. The prac¬ 
tical approach is the one described above, namely, to make the basis successively narrower. 
In that process, the number of bound states will increase, and the states apt for an EFT(^) 
description are those which are, one, relatively close to threshold, and two, stable under the 
addition of narrower basis states, f.e., independent of Arqm- All graphs show the lowest 
breakup threshold set by B{a). At NLO, this threshold is at the experimental a binding 
energy independent of A as it is included in T)'. For A = 300,400 MeV, predictions were 
made with interactions whose three-body parameter was fitted to B{a,exp) at LO. Hence, 
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FIG. 18: Predictions for the shallowest six-nucleon energy eigenvalue in the six-nucleon system with 
EFT(]t) as a function of A. The (blue) gray band shows (N)LO uncertainties assessed via cutoff and 
input variation. 
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the threshold is at the same position for (N)LO and for all three A values. As the potential 
with A = 200 MeV uses B{t, exp) to hx A 3 N 1 at LO, B{a) (gray line) is a prediction at this 
order. 

For each interaction, we tentatively identify the bound eigenstate closest to the thresh¬ 
old with the EFT(y') prediction for the ®He system. At LO, bound states (empty 


squares in hg. 18) with the following energies are predicted with the 720-dimensional 
basis: i?Lo( 6 ,A = 400 MeV) Fa 55 MeV, i?Lo( 6 ,A = 300 MeV) fs 35 MeV, and 
Bi^o{6,A = 200 MeV) fs 43 MeV. This implies a LO uncertainty of about 20 MeV 


(gray band in fig. 18). At NLO, the corrections place the states (filled squares in 


fig. 18) at: i?NLo( 6 ,A = 400 MeV) fs 47 MeV, i?NLo( 6 , A = 300 MeV) fs 41 MeV, and 
i?NLo( 6 , A = 200 MeV) fs 44 MeV. The pNLO contribution to the binding energy is attrac¬ 
tive for the smaller two cutoff values and repulsive for A = 400 MeV. EFT(^) predicts thus 


at NLO a binding energy for the state between 47 MeV and 41 MeV (blue band in fig. 18). 
Therefore, for this state, the theory does not converge to the physical, shallow ®He state. 
Yet, this result does not necessarily mandate a promotion of a four-, five-, or six-nucleon 
interaction for an EFT(^) description of the six-nucleon system which converges to exper¬ 
iment. If such an unusual promotion of a momentum-dependent six-body vertex would be 
necessary, a behavior of the six-nucleon spectrum similar to that of the three-nucleon spec¬ 
trum is expected following a change in the renormalization scheme without the additional 
constraint that fixes, e.g., B{t). In the six-nucleon sector, the missing constraint should 
result analogously in a shift of the spectrum (limit cycle behavior). Instead, the considered 
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eigenvalues are found within defined uncertainty ranges which decrease from LO to NLO. 
This observation suggests that the RGM basis was chosen wisely for the states investigated 
here, but poorly for the physical states. 

First, it could be insufficient for an expansion of a shallow state bound by a relatively 
small amount below threshold. Specifically, a much broader range of width parameters 
for the expansion of the two radial functions of the halo-neutron coordinates ni ^2 has to be 
considered in order to discriminate such states from a-nn continuum states. Misinterpreting 
such a free-a state as a bound six-nucleon state amounts to a four-body calculation without 
predictive power. To illustrate this point, we included the predictions for such a state in 


fig. 17 for A = 400 MeV (filled and empty circles, bottom right graph). An unbound a-nn 
state of about 24 MeV is there not to be confused with a bound ®He at NLO, although 
the pNLO correction places the state very close to threshold. If the width parameters for 
ni ,2 cover a range beyond the expected halo extension, the numerical ambiguity between 
a-nn continuum and ®He bound states can be resolved as follows. The continuum states are 
naturally much more sensitive to changes of the basis regarding vectors that correspond to 
very large separations of a neutron from the core. The problem is to induce such changes 
without affecting the expansion of an also widely extended shallow bound state. 

An example of such a modification in the large model space is shown in the upper graph 
for A = 300 MeV in fig. Two bases, each of dimension 2436, were chosen which differ 
only in the width set used to expand the dependency on the coordinate of the second halo 
neutron. In one basis, the smallest width of this set was replaced by an order-of-magnitude 
larger value. In effect, this more localized basis expands the ground state (squares) with 
B{6) > 60 MeV as well as the broader basis but deviates significantly in its prediction of 
the shallower state (compare triangles (“more localized” basis) to circles (“broader” basis)). 
A similar sensitivity to the removal of the broadest basis vectors will eventually only be 
observed for the a-nn continuum states if those vectors are sufficiently broad to have neg¬ 
ligible overlap with the halo states. With the resources available, this investigation cannot 
be included in this work. Restating the problem, in general one deals with the slower con¬ 
vergence regarding the basis size of a variational calculation for excited states. We propose 
here that the ®He halo state is an excited state in the spectrum of the employed poten¬ 
tials, in contrast to the deuteron and triton. For the latter two systems, the regularization 
through the EFT cutoff A or the model space Argm did not allow for (admissible) deeper 
lying states. In the six-nucleon sector, the chosen regularization parameters admit such 
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states and complicate the numerical analysis. The three-body analog would be to identify 
the triton as an excited Ehmov state for an interaction which allows for more than one 
bound state. The other states in this spectrum are unphysical. Yet, the formal approach 
is the same as the one where those deeper states are observable, for instance in the 
isotope. For this nucleus, the possibility of an excited ^®C-2n halo state was investigated!^ 
in an effective three-body framework. 

In contrast to this hypothetical carbon halo, the P-wave resonances in the a-n subsys¬ 
tem are supposed to dominate the 6 -helium ground state. In the previous chapter, no such 


resonances were found at LO (see fig. 13). A future scattering calculation will determine 
whether only interactions which do predict an a-n P-wave resonance also result in a shal¬ 
lowly bound 6 -helium. The results presented here would be consistent with this finding. 
Here, we employed only potentials without such a hve-body resonance. As a consequence, 
the only bound six-body structures found are unphysical, i.e., the system is unbound rela¬ 
tive to the a-nn threshold. The EFT does not claim to predict whether the state is bound 
or unbound. Assuming uncertainties for the shallow state of the same order as assessed 
here for the deeper states, even at NLO, higher-order modifications (of order 7 MeV, see 
above) could push the state above threshold. 

Before I conclude, the limitation to small cutoff values relative to the ones used in 
A < 6 calculations, and how it can be overcome in future analyses is explained. For a 
simplification of the numerical treatment of larger cutoff values, I propose to increase A to 
values for which the a is no longer the ground state of the resultant potential but rather 
the shallowest state of a tower of four-body bound states. This would utilize a phenomenon 
observed in the three-body sector where “new” bound states enter the spectrum for discrete 
values of A at threshold. Increasing A from such a value also renders the initially shallow 
state steadily deeper and deeper, until the next bound structure enters, again at threshold 
(see e.ff. ch. 3.4 in HI)- Analogously, I expect additional four-body bound states entering 
the spectrum - either from above or below - at the ^He-n threshold at certain cutoff values 
At. Increasing A from such a At until the lowest state is bound by approximately B(a, exp) 
and identifying this instead of the, by then, much more deeply bound state(s) with the a 
would also mandate only smaller pNLO corrections. One would monitor B(t) and B(a) 
simultaneously for A 3 N 1 = 0 and select cutoff values for which the 3NI has to contribute 
relatively little in the course of renormalization. This avoids cutoff ranges where the LO 
renormalization requires large modihcations from the 3NI. The complications associated 
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with large 3NI LO modifications can be understood by considering the demands on the 
variational basis. Two states of very different spatial extent have to be expanded: the 
unrenormalized deeply bound and thereby localized state, and the more extended states 
renormalized to the physical B{t) and B{a). Initially, I saw no need to draw this option, 
and used A values as given above. 

In conclusion, the convergence of relatively deep states in the six-nucleon spectrum from 
LO to NLO strongly supports the assertion that a shallow state would be predicted by 
EFT(it) without a promotion of a higher-order vertex. The finding is in contrast to the 
limit-cycle behavior in the triton channel of three-nucleon system. The identification of the 
discovered states with ®He would be false but their existence and convergence are hints of 
the existence of a more shallow state in the spectrum. Whether this state was not resolved 
due to the small RGM basis, or the fact that there is no bound 6-helium for the employed 
interactions, but rather a shallow resonant state, needs to be addressed in a future larger- 
scale bound and scattering computation. 

X. SUMMARY 

Nihil est dictum, quod non sit dictum prius - Advances in the understanding of particular 
areas in the parameter space of the physical world - to recapitulate this object from ch. |I] 
where the effective field theory formalism was introduced - are guided and inspired by 
experience in other areas. This process of generalizing theories - developed to describe a 
particular system - to others, refining particular aspects while discarding other, irrelevant 
components of the theory, is more efficient if the principal path between the areas is known. 
Technical advancement makes larger areas of the parameter space available for scientific 
investigation, and by that discovers paths connecting one area of nature to another where 
the understanding might be either more mature or still in relative infancy. In nuclear 
physics, this general perspective focuses on three current areas of interest. 

The first is the relation of the nuclear interaction to its underlying theory of QCD. 
Recent developments in calculating nuclear properties from this high-energy theory include 
the prediction of the neutron-proton scattering lengthen, the nucleon masses!^, electric 
and magnetic moments of nucleons® and even tri-nucleon binding energies®. 

Similarly, taking the NN interaction as fundamental, its systematic extension to describe 
many-body systems in terms of a density functional theory® falls into the category of 


85 


relating different parameter subspaces of the physical world. Both developments depend on, 
but also drive, the advance in EFT methodology (xPT, ^EFT^Sl, EFT(y')) and numerical 
methods (lattice techniques, quantum Monte Carlo). 

While in these areas increasing computational power makes certain systems accessible 
for theoretical investigation, experimental skills in controlling atomic systems at very low 
temperatures allow the investigation of universal properties of bosonic® and fermionic^i^ 
systems. The modification of Feshbach resonances® of trapped atoms via magnetic fields 
results in a change in the atom-atom scattering length uaa and measurements of loss rates 
out of the trap yield information about dimer, trimer, quadrumer, etc., formation as a 
function of qaa- By that, aspects of atomic physics are related to nuclear physics, where 
similar experimental techniques are nonexistent. 

This sets the context for this work, where I have presented a study of an effective field 
theory without pions in few-nucleon [A < 6) systems. The EFT formalism was used to 


construct the interaction (ch. Ill) and to refine it model-independently. The EFT(]t) is the 
low-energy derivative of a pionful theory and thus linked naturally to QCD. 

We demonstrated explicitly the usefulness of the theory by presenting results for binding 
energies of the triton and a particle, the charge radius of the triton, and the S-wave neutron- 
deuteron scattering length in the triton channel of increasing accuracy from leading to next- 
to-leading order. The first-time (except for the nd scattering length and B{t)) application 
of NLO EFT(^) to those observables conhrmed assertions based previously for ^ > 4 on LO 
calculations only, which I summarize: the 4-helium bound state is a universal consequence of 
and neither for the three- nor for the four-nucleon system, the power counting had 

to be modified by, e.g., the promotion of a four-nucleon interaction to NLO, in order to be 
consistent with a priori accuracy predictions of EFT(^). More specifically, the effective field 
theory formalism can be used to explain empirically found correlations amongst few-nucleon 
observables like the Phillips and Tjon line. In these two cases, it relates the deviations from 
data to, first, an incomplete renormalization in the three-nucleon sector, and second, to 
higher-order interactions omitted in the course of the EFT expansion. The theoretical 
uncertainty at every order of the calculation can be quantified in this approach. Numerical 
inaccuracies were demonstrated to be negligible (see discussion of figs. &: [T^ . Different 

results for different short-distance parameterizations therefore have other origins. The size 
of the expansion parameter Ptyp/Abreakdown determines how fast an EFT expansion converges 
and decides its usefulness for the calculation of an observable in a given system. In heavier 
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systems like 4-helium, the expansion parameter can a priori be as large as 1, so they are 
border line. However, we find that next-to-leading order corrections to leading-order results 
are still parametrically small in the four-nucleon system. This confirms a pattern already 
seen in two- and three-nucleon systems at momenta which approach the a priori breakdown 
scale but where convergence is still found, see e.g. [2^ 1^ 1^1^ . 

The feasibility of calculating few-body observables with the numerical Resonating Group 
Method and EFT(^) potentials was thus established. The RGM has so far been employed 
in conjunction with potential models of the nuclear {e.g. the AV18+UIX/IL2 predictions in 
this work for ao(^He-n)) and quarkEH interaction only. Rs practicality to obtain predictions 
with the very different radial structure of the EFT(]if) interaction - single Gaufiian functions 


(eq. (26)) - was not guaranteed. The reliability of the method was assessed in a neutron- 
deuteron scattering calculation, where the predictions of the RGM agreed with the Faddeev 
results and also showed the impact of different input data on the estimate of LO uncertainty. 
The conclusion from this analysis, in combination with the correlations of bound-state 
properties for A = 3,4 (fig. 10) is that for an estimate of the renormalization- 

scheme dependence of observables the data input and the regulator should be varied. 

For two correlations for which leading-order calculations exist, namely between the bind¬ 
ing energy and charge radius of the triton (fig. and the binding energies of 4-helium and 
the triton (fig. [To|), our coordinate space EFT(/) calculations at NLO report the expected 
improvement from LO to NLO consistent with an expansion parameter Ptyp/-^breakdown ~ |- 
By that we demonstrated that a consistent description of the a particle is possible at NLO 
in EFT(if). In fig. we also report a correlation between the triton binding energy and its 
difference to the ^He binding energy. As the EFT(^) potential is at NLO isospin-symmetric, 
this model-independent difference is attributed entirely to Coulomb interactions, which are 
included in the RGM. At the physical triton binding energy, this value agrees well both 
in magnitude and uncertainty with estimates of charge-symmetry breaking and Coulomb 
contributions to ^He binding. In fig. 


II 


a new correlation between the triton binding energy 


and the real part of the singlet S-wave scattering length of 3-helium-neutron scattering sim¬ 
ilar to the Tjon line is presented. This, and the three aforementioned correlation bands, let 
us also conclude that no four-body contact interaction is required to renormalize the system 
at next-to-leading order. The position of all four bands, which represent universal proper¬ 
ties of the two-nucleon system, was determined by fitting nine NN potentials differing at 
short distances but with identical long-distance behavior, by variations of the 3NI strengths. 
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and by changing the cutoff. Consistent with a basic tenet of EFT, namely model indepen¬ 
dence, we also showed that the results of the phenomenological models AV18(+UIX/IL2), 
which share the input of our EFT(^) NLO potentials, agree with their results within NLO 
accuracy. 

The first-time exploratory application of EFT(y') at LO and NLO (perturbative) to the 
5- and 6-helium systems considers low-energy a-n scattering (LO only) for a large range of 
cutoffs, and bound-state calculations for both systems at LO and NLO. The renormalization- 
scheme dependence for LO EFT(y-) for the a-n S-wave scattering phase shift was assessed 
thoroughly and was found not to exceed the naive 30% deviation from data even for A —)■ oo. 
The accurate prediction of i?(a, exp), either by an appropriate choice for the cutoff or by 
enforcing it as a condition when fitting the LECs, did not reduce the uncertainty. 

The five-nucleon system was found unbound in the ^He channel at LO and pNLO. A 
calculation without the Coulomb force led to the same conclusion that an unbound 5-helium 
system is a universal property of the nuclear force. 

The ROM bound-state LO and pNLO EFT(^) calculations in the six-nucleon system 
identified deeply bound structures. The uncertainty in the binding energy was found to 
decrease from LO to NLO, but no agreement with experiment within the assessed NLO 
uncertainty was observed. The absence of a shallow bound state in the spectrum and the 
convergence of deeper lying states was interpreted as either a shortcoming of the numerical 
method or an unbound shallow resonant state. The latter is deemed to be the more likely 
scenario, because we could not reproduce P-wave resonances in the ^He subsystem. The 
results do not contradict EFT tenets. The observed decrease of the uncertainty in the six- 
nucleon binding energies rather supports the assertion that no higher-order operator is to 
be promoted for an EFT(^) description of a shallow six-nucleon state. 

We argued that a larger RGM basis is imperative to test the applicability of EFT(^) 
to the shallow states in both, 5- and 6-helium. By considering NLO shifts which leave the 
state bound, this can be done in the computationally less expensive bound-state framework, 
while the general case of a transition from a shallow real- to a virtual bound state and vice 
versa, can only be analyzed in a larger-scale scattering calculation. This work also provides 
the feasibility study for this forthcoming investigation. 



Appendix A: Next-to-leading Order Potential in Coordinate Space 


The transformation from a potential in momentum representation {p'\V\p) to coordinate 
representation, which can be used in our RGM implementation, is given in this section. 
Without specifying the form of the potential operator, the form of interest is given by: 


(r-'lPlT) = 


(27r)3 


d/r (r |l/|r)(r|T) 


(Al) 


Changing the basis yields: 

1 


{r\V\f) = 


d^p'd^p {r'\p'){p'\V\p){p\r) = 


(27r)6 
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1 


(2vr) 


d^p'd^p e^^-^Vi^,p)e-"-P 


d?kd^q 

d^kd^q . 


(A2) 


The operators in V^^^^'^{q,k) (eq. (15)) are polynomials in k. Replacing k = in the 
potential operator with the derivative acting on the exponential only, yields an equivalent 
equation. This replacement must not be confused with replacing k by its coordinate space 
representation, in which case, r would be the conjugate variable to k. 

Consider, for example, the operator multiplying in eq. (15), 

= e ^ di ■ k d 2 ■ k. 


{f^\V\f) = J d^q e^^<^^W{q,iVr) j d^ke^^<^-^ 

= - ai^ia2jVr,iVr,jS{f' - f) . (A3) 
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This operator acts on a wave function, T(r) = 


(f'|T|T) = -cti,jCT2j / d^r/o 
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(72 • rdi ■ V 


(A4) 


^ 1 A^ _ _ _ \ 

+CTi • r(J2 ■ V + -CTi • (72- —di ■ 7(72 • 7j 4/(r) 

Derivatives are understood to act on the expression in the first square bracket on their 


right, i.e., on the deita function in eq. (A4) and not on the wave function, because they 
were introduced to ’bring down’ the /c’s from the exponentiai. Anaiogous transformations of 


aii operators in eq. (15) aiiow for an expression of the potentiai in coordinate space as given 


in tabie XI Hermiticity is expiicit in this form, and the anticommutator structures aiiow 
an efficient anaiytic evaiuation of the corresponding matrix eiements reiative to a form iike 


nai’veiy derived in eq. (A5). 
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Appendix B: Genetic Search Algorithm 


With the function chosen to regulate the contact interactions (eq. (20)), an analytic 
calculation of the LECs, even at leading order, is not possible. Furthermore, a certain A only 
suggest how select and does not provide a rigorous way to derive values for the boundaries 
of an interval from which the Gaufiian width parameters 7 „, which define the RGM basis, 
are to be taken. Optimized values for both the LECs and the 7 n,’s are determined here with 
a genetic algorithm ES . This stochastic method was chosen because: 

• The entire parameter space is probed. This made the determination of nNLO LECs 
feasible even without an estimate of their magnitude. Such an estimate, available for 
the LO LECs through eq. (13), would be based on the naturalness of the 


only, and therefore be unreliable because the LECs are redefined. The Aj’s in eq. (26) 
are linear combinations, Ai = hence a change in the original 

C'j^N)LO’g percent level could induce relatively large changes in the Aj’s because 

of the A dependence. 


• The search parameters can be adjusted to reduce the probability of identifying a 
local instead of a global minimum. Local minima could, for instance, correspond to 
an interaction with a deep two-body bound state in addition to the relatively shallow 
deuteron. This ground state, however, is beyond the range of applicability of the 
EFT, i.e., Argm might not even allow for its expansion. In general, its energy will 
depend nonperturbatively on all renormalization parameters. 


• The computation factorizes into independent blocks and scales almost linearly with 
the number of processors used. 


The task can be abstracted to the problem of finding a set of parameters {a,} that minimizes 
the function W ({ui}). Specifically, 


class (i) In a given complete model space, the LECs are to be fitted to data, i.e.. 


W({a,]) = E '"‘TTxpT**' (“■) = 

(implemented for this EFT(^) application). 


{Ai,,,., 9 } and V as defined in eq. (28) 


class (ii) For a given potential, dimension, and angular momentum structure, the 7,1 
are optmized to maximize the binding energy, i.e., VE ({oj}) = and 

{ai} = { 7 n} (implementation from [99]). 
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The incarnation of the algorithm used here proceeds along 8 basic steps: 


1. A set is mapped onto a bit representation which is called individuum. 


Oi 


(Xi 






On 


i|i|i|o|i|o|o|o|o|"ir 


In this work in particular, Oj’s were either LECs A, or RGM basis width parameters 

7- 


2. A random set (initial population) of Z individuums is generated from a set of bound¬ 
ary conditions of allowed values. Every individuum is assigned a fitness value through 


the function W {e.g., eq. 28). The individuums are ordered according to this fitness 
value Wi 


Cl 




Cn I 


^1 I I I I I 


ai 




On I 


'^1 — '^max 


strongest individuum, 

highest probability of becoming a parent; 


Wz — 'd^min weakest individuum, smallest probability to procreate; 


with Wi < Wj Vf < j 


3. Two individuums (parents) are selected from the population. The rank (6) of the two 
parents in the hierarchy introduced in the previous step, is determined by a random 
number x G [0,1] and b = x-\J^ln{x) • f ■ If 5 < Z, the element b is one chosen 
parent. This procedure ensures a higher probability of selecting parents with a small 
W value but does not entirely neglect entries with high W. 


4. The two parents are used to create two individuums (offspring) via uniform crossover 
(equations for bits Xi and Zi) with a random individuum m (mask). 



mask 




Xi = {wi A mi) V {ui A mi) 

Zi = {wi A rhi) V {vi A mi) 


offspring 


Uniform crossover is only one of infinitely many choices for a rule to create offspring. 
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5. Each offspring is mutated, i.e. every bit is inverted with a probability of p. Here, 
p = 0.01 was used. 

6 . The W values of the offspring are calculated. 

7. The offspring is inserted in the population according to its W value. 

8 . The two individuums with the highest W value are removed (only the strong survive) 
from the population. 

Steps (3) to (8) are repeated k times. This iteration together with the preparation phase 
(1), (2) is called run. At the end of a run the individuum with the lowest W value is 
returned. Consecutive runs were used, each time only the intervals - set in step (2) - were 
adjusted. At the end of a run, they are centered around the values of the best individuum 
of that run. Additionally, the size of the intervals is reduced by a factor of |. After that, 
the next run is started with the new bounding intervals. Typical populations consisted of 
Z ~ O individuums, which created their respective Vitruvian Man within less than 

O (lO^) generations. 


Appendix C: How an Integral is done 


The analytic expressions for a coordinate space matrix element of the operator oc Ag in 


y(NLO) (26)), in the bound-state basis (eq. (|35[)), i.e., set parameters {f,j,lj,Sj,jj} 
in bra and ket, are derived in this section. The operator structure is given by: 




r,T' uiy' 


with 


(Cl) 


Vjj(^) ly : spherical component U of the derivative 
acting on the relative coordinate between 
particles Z and J Jacobi coordinate T' , 

PV,, = , A(') = . 


(C2) 
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The integral J{V; Li, Li^,... ,Lz,Lz^) (eq. (45)) is obtained from a generating integral 


which is defined in terms of the generating function of the spherical harmonics (eq. (50)) 


by: 


...,a„b,) = J d^si ... d^SA-ie . 


~ y~] / ^nbn'Qn 

,/g MM ^ ^ ^ 


E n 

\n=l 


7 I 

bn. • 


J'^{'P-,li,mi,... ,lz,mz) (C3) 


The number of particles in the ket state is ric^. For all states used in this work, this number 
is identical to the corresponding number rici ™ ^be bra, and defined here only as a reminder 
that a grouping of nucleons in S-wave clusters within a fragment, leading to possibly different 
coordinate structures in bra and ket, is an unused option. The exponential of the interaction 
operator, exp( —^f^) is implicit in the transformation matrix for each summand 

of the anticommutator. A principal axis transformation. 


A-l A-1 A-1 A-1 

and Qn = Pnfi'tf,' such that • V = X] ’ 

A=/i p'=l A 

(C4) 


followed by a translation to complete the square in the exponent, 

y^.r, CLnPnXbn 


t\=Ux + 


2/3a 


(C5) 


results in an expression of the integral in terms of the norm integral, defined analogously 


to the operator generating integral in eq. (50): 


I ,h\, ... , ttz ,bz) — .f^Norm (®1) ) • • • ) ) bz) ^ ^ ^nn' ■ CLnO-n' Pq (^b) 


with 


7 _ .,AA' 

7X71^ / ' Tt' 


A-1 


{Pr'fi + P'Iit') {P%' + p'li'r) 


L 

A-1 


E 


“ PnaPn'a' + 0{nc, - 1 - n)9{nc,. - 1 - n')?n'r?nr' 


1-1 

E (Pri^ + - 1 - ACnr' 


a,a 
A-1 


A-1 


11=1 


a=fi 


TfiaPn' a 

Pa 
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and 


P,= < 


1lu2 
3 


2Z)2 Z,2 


q = 2 

2{bnbn' + bl,bn) 


q = 1 

1 + ^bnbn' + b^i) 

for 

z? = 0 

— 2{bn + bn') 


g = -1 

2 


q = -2 


Relating the coefficients of anh^ and anbl^a^,b^', in eq. (C3) and eq. (C6) gives the matrix 


elements for the five spherical components of the operator: 

/ / 

jij,q=2 _ \ '' o ^ TL n' ^ ^ jNorm 


I ^ ^ 2,771^1 + 2 yNorm 

“T 2 _^^^nn 




a,7 




n^n' 




^ rNorm 

'-^Zn-l,mn+l'-'Z„,-l,m„,a(P;..,,z„_Pm„+l,...,Z„,-l,m„,,..) 


rNorm 

TNorm. 


(^n l)^nC'Zn-2,m„+l J-Norni 


^ln,r. 


jihq=0 _ ST' ly 

Y 3^™' 

n^n' 




/~i rNorm 

rNorm 

+ ^(6 ~ jNcm 





(T’;...,Z„- 2 ,m„....) 7 


Tij,q=-^ 

{l,m} 


^(-2)Zn 

n^n' 


Ifiin' 


rNorm 


^ rNorm 

OZ„-l,m„<--Z„,-l,m„/-l-'(P;,..,Z„_l,m„,...,Z„,-l,m„,-l,...) 


+ ^(-4)Z, 


{In l)^nC*Z„-2,m„-l j-Norm 




('P;...,Z„-2,m„-l,...) ’ 


Tij,q =-'2 

{l,m} 


n^n' 

z 

+ 2Z„, 


Inin' 


{In l)^nC^/n—2,mn—2 


C /O rNorm 

Z„-l,m„-lc^Z„/-l,m„/-l«^( 


{'P;...,Z„-l,m„-l,...,Z„/-l,m„/-l,...) 


c, 


In.Tnn 


rNorm 

'^(P;...,/n-2,mn-2,...) ’ 
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with Clm defined in eq. ( [^ and {l,m} = ... ,lz,mz)- Given the values for the 

norm matrix elements®^, the operator elements are thus obtained for permutations V 
representing double cosets®], only. 


Appendix D: The pNLO calculation with the RGM 


Thissection documents the modifications to the adopted RGM implementation in order 
to carry out perturbative NLO calculations as explained in ch. H The stage of the RGM 
implementation at which the extensions were made and the technical notation used here is 
reported in [49] for bound-state calculations and in unu for the scattering problem. 

Obtaining perturbative corrections to bound state energies was already implemented in 
the program chain, while the calculation of the operators with coordinate space structure 
and {e~^^ is new. The operator was implemented as the fourth operator in 

the LUDW part of the program chain. For this operator, LUDW was modified to prepare the 
calculation of reduced spacial matrix elements (see eq. @) for 

=-\/3-— — , (Dl) 

3 dvr ^ 

'-V-" 

=F 

with spherical components of the relative coordinate between interacting particles equal 
to solid spherical harmonics®], r'^ = \J~^ ■ yiq{r), and a prefactor F implemented in the 
KOBER part of the program chain. 

The peculiar anticommutator structure chosen to represent the other momentum- 
dependent NLO contribution (compare to the form in [102] i made it possible to calculate its 
matrix elements by a linear combination of modified matrix elementsof the kinetic energy 
operator. The changes were made in the QUAF code, in which the width parameters of the 
basis and the interaction are specified. The derivation of an analytic formula for a generic 
matrix element parallels the one given in app. 0 In comparison to the more complicated 
case of an operator of rank 2 as considered in app. the structure for pNLO is a scalar: 






(D2) 


The two summands of the anticommutator yield expressions for their matrix elements in 
the RGM basis which differ only in the definition of the transformation matrices: o' , is 


^-2 ^ ^ 

modified for the e element and p"^^, for (see eq. (C3)). This ingenious 




97 









treatment [104) makes the explicit calculation of operators of the form r • V and • fa ■ V 
unnecessary. 

With those extensions, all numerical tools are available to realize the NLO constraints 


(iii)-(v) as defined in ch. Ill, In distinction to the ground state wave functions for the 
deuteron, triton, and ^He and the asymptotic fragment relative functions, the scattering 
states (r|'I'^(fe)) are complex valued. The reality of a bound-state wave function (r|d,t,a) 
is a consequence of the real RGM basis and the real potential, while a boundary condition 
enforces real valued relative motion radial wave functions between the fragments. In the 
hf basis (eq. (36)), the NLO amplitude - given diagrammatically in fig. [^- is 


{h-\f\h+) = {h-\no\ht) + 


(D3) 


with 


T + Rlo 


l^h) =- 




The wave functions (riT)}") are obtained by calculating the S-matrix in LO ^lo (see steps 


leading to eq. (38)). Substituting the Gaufiian expansions for hf in eq. (36) and the LO 


are 


S-matrix then yields the expansion coefficients for (riT)}"). The coefficients for (r|'I'Q 
calculated by (rlT))") = A modified version of the SPOLE component of the 


program chain writes an input file in the standard shown in table XII This file specifies: 
the functions ) and Gi,Fi in terms of complex expansion coefficients; S'lo; the energies 
for which the states were calculated. The information is processed by an updated DR2END 


code which writes the pNLO S-matrix to an output file with standard as in table XIII The 
LO and NLO part of the S-matrix, 


-Snlo = 'S'lo + 




V 2 ^ 


) (cce- 


9 A- 

CrV e 4 


-I- Cp{e 


\vn)- 

(D4) 


are printed separately. 


For LECs Cc,r,p of order 1, the 5 nlo matrix elements between the scattering states, 
with expansion coefficients Ci('I'^) determined in SPOLE, typically assume values between 
10“^ and 10“®. Gomparing magnitudes of real- and imaginary parts of the S'nlo matrix 


elements (eq. (D4)) for: the three NLO two-body operators, for a single operator, and for 


individual operators the matrix elements in the ^5*0- and channel, differences up to 
10^ are common. The discrepancies are increased by the regularization of Gi (eq. (39)) 
which reduces the size of the matrix elements. Increasing the parameter (3 - defined in 
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Ne 




Re{fcl = -v/2TOredAl} 

Imjfci} 



Re{/i(A:i)} 

Im{/i(fci)} 

Re{ffi(fci)} 

lm{gi{ki)} 





Re{/v(fci)} 

Im{/Ar(fci)} 

Re{5rAr(fci)} 

Im{gAr(fci)} 

Re{*S'Lo(5i)} 

Im{5'Lo(5i)} 



Re{ptiki)} 

Im{p+(fci)} 







Re{p)^(fci)} 

Im{p+(fci)} 



Re{pr(^i)} 

Im{pj-(/ci)} 







Re{Pw(fci)} 

Im{p^(fci)} 




TABLE XII: Format of the modified INEN file generated by S-POLE_DWBA.f and processed 
by DR2END_DWBA. f. The head of the file remains unchanged with NBAND1=10(13) for a 
NN(NN+3NI) calculation and NZZ=-2 to read the LO data. The additions are attached to the 
end of the file. The coefficients expand the functions: GAkr) = \ , 

= \J and (r|^',=‘=(fc)) = (*)e"^'’'''with width pa¬ 

rameters 7 i set in INQUA_N. 


Reffcl = \/2>71redR'l} 

Re{S(fci)} 

Im{S(fci)} 

Re{5Lo(fcl)} 

Im{SLo(Ai)} 

Rel^NLof^i)} 

Im{SNLo(fci)} 








Re{fc]v^} 

Re{S(fc^J} 

Im{S(A:iVe)} 

Re{5Lo(fciVe)} 

Im{5Lo(A]Ve)} 

Re{5NLo(fc]Ve)} 

Im{5NLo(fcjVe)} 


TABLE XIII: Format of the output file generated by DR2END_DWBA.f. 


eq. (39) - to allow for a non-zero Gi at smaller distances, and therefore a larger overlap of 
the suppots of the wave function and the Gaufiian-shape potentials, comes at the expense 
of a less accurate expansion of Gi in the Gaufiian basis. Instead of increasing the basis 


size, I chose smaller cutoff values A for which the S'nlo integrals (eq. (D4)) with expanded 
functions Gi were practically identical to the one calculated with the unregularized Gi. The 
parameter /3 constitutes another part of Argm and was set to (5 = 1.1 fm“^. Although, the 
insensitivity of NN phase shifts to this parameter has been shown, the results found here 
for five- and six-body systems inferred from interactions based on the newly implemented 
pNLO calculation, require for a new investigation. In addition to the suggested refinement 


of the RGM bases in ch. VIII &: IX, the /3 dependence of 5(5,6) has to be assessed. In 
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other words, insensitivity towards high-energy physics parameterized by /3 remains to be 
proved. A priori, there is no reason why /? should be less significant in this respect than 
the truncation of the basis by a finite set of width parameters. 


Appendix E: Ansatz for the triton wave function 


In this section, we use the triton system as an explicit example for the construction of 
the RGM trial wave function. The chapter is thus an amendment to the general discussion 


in ch. IV following eq. (42). 

The proton and the two neutrons can assemble into three different two-fragment groupings: 
deuteron-proton {a = 1), singlet-deuteron(d)-neutron (a = 2), and (di-neutron)-proton 
{a = 3). For a = 1,2, the definition of the Jacobi coordinates pi ^2 is identical but the 
structures differ in the intermediate coupling 512 of the spins in the two-particle fragment 
(see fig. 19). Using the notation introduced in eq. (42), the explicit ansatz for the antisym¬ 
metrized (^) variational wave function with total angular momentum J = ^ reads: 


1 

)2 = 


Ay.' 




s Q=1 



(El) 


The square brackets indicate an angular momentum coupling as defined in eq. (41). The 
spin functions H corresponding to the respective groupings are given as: 


‘=‘1 ~ 


‘=‘2 ~ 


‘=‘3 ~ 


[np] n 


[np]*^ n 


r lO 
nn n 


(snmiSpm2|lm3) (lm3S„m4|l/2m) |s„,mi)|sp,m2)|sn,m4) , 

™1,2,3,4 

(E2) 


im 


where the individual nucleon spins \sn,p = \,m) are coupled to a total spin 5 = ^ equal 
to J because the orbital angular momentum is zero. The product states are superimposed 
with standard Clebsch-Gordan coefficients. The ansatz for the coordinate-space component 
of the wave function, assuming S-waves only, is: 


^ Caa'^aa' = ^ C^a' exp { -71 (aa')/?? “ '^2{aa')pl] ' [3^0 {Pl) ® 3^0 {pAif ° , 

OlOl' OlOl' 

(E3) 

with a total orbital angular momentum I = 0 and variational parameters Caa'- A basis 
vector is then labeled by its grouping a and a specific pair of width parameters 717 indexed 
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(a = 3) 



FIG. 19: Three-cluster fragmentations of the triton. The left panel shows the dominant grouping 
including a bound deuteron fragment with intermediate spin S 12 = 1 and the singlet deuteron with 
512 = 0. The structure on the right shows the nn grouping which allows for 5i2 = 0 only. 


by a'. The convention of |l 8 ] is used for the solid harmonics y. With the trial wave function 
in this form, as a linear combination of nonorthogonal basis functions, 


/ = [^c 


P - 
0 . 0 . 


(E4) 


the (Ritz) variational method reduces to an eigenvalue problem of the Hamiltonian inside the 


space spanned by the (l)aa'- In table XIV a toy basis is defined with four a = 1, three a = 2, 
and two a = 3 vectors. The widths were chosen of similar size to represent either compact 
or extended distributions. For example, the a = 1 set 71 = 11.1 fm“^, 72 = 8.2 fm“^ has a 
very small support relative to the vector with 71 = 0.9 fm“^, 72 = 0.1 fm“^. The former is 
needed to refine the approximation of the wave function at short distances, while the latter 
would be more important to expand extended structures or scattering functions. 

With a LO interaction of the form 

3 


H (-40.6 - 4.8 <Ti ■ aj) . 


(E5) 


i<j 


and a variational basis of dimension d = 1 to d = 9, the generalized eigenvalue problem, 

d 


^((4i^i#)-ev(3)(4i^I0. 

i=l 


Cj = 0 


(E 6 ) 


has to be solved. To simplify the notation the pair aa' is replaced by a single index. 
From this equation, the lowest eigenvalues as listed in the fourth column of table |XIV 


are 


calculated. The basis is defined by eq. (E3), and widths are given in table XIV where d = 1 


uses the width pair of the first line, d = 2 the first and second line, etc.. Adding vectors 
to the basis which correspond to highly localized structures relative to the range of the 
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TABLE XIV: Width parameters defining the variational space for a toy triton when used in eq. (E31. 
In the fourth column, the displayed smallest eigenvalue of the Hamiltonian was obtained with a 


potential of the form eq. (E5) and a basis which includes all vectors above the respective line. 



a 

a' 

7i [fm ^1 

72 [fm 

EV(3) [MeV] 

d-n 

1 

1 

5.9 

3.5 

480 


1 

2 

0.9 

0.1 

10.5 


1 

3 

11.1 

8.2 

10.5 


1 

4 

3.3 

0.5 

9.48 

d-n 

2 

1 

0.7 

0.4 

4.48 


2 

2 

1.1 

0.4 

3.34 


2 

3 

16.1 

11.1 

3.33 

nn-p 

3 

1 

3.3 

4.2 

3.33 


3 

2 

0.9 

0.28 

2.48 


TABLE XV: Expansion coefficients Caa' for the three-nucleon ground state for a basis of dimension 

9 as defined in table llOVl Row number labels a and the column a'. 

12 3 4 

110.65223 0.51435 -1.7785 -0.54223 

2 -1.7899 1.6454 -1.0757 

3 0.63488 0.47931 


potential - line 3 and 7 in the table - have no significant effect on the eigenvalue even spaces 
of such a low dimension. In d = 9, the eigenstate with the lowest eigenvalue EV(3), as given 


in table XIV, is, for example, expanded with coefficients Cj (eq. (E6)) as given in table XV 


As fragment fl in a triton-proton scattering wave function, this state would be in the 


channel wave function as given in eq. (40). 
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